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ABSTRACT 


In order for railgun technology to be relevant to the Infantry, the design of the 
projectile must incorporate the following three concepts: an effective ballistics package, 
geometries for aerodynamic stability, and a non-parasitic conducting armature. I 
designed an effective 30mm and scaled 40mm projectile that incorporates the 
aforementioned concepts. My ballistics analysis concluded with two AUTODYN™ 
finite-element computer models that refined theoretical estimates for target penetration. 
The proposed railgun projectiles were effective in penetrating 100 mm of Rolled 
Homogenous Armor and in perforating 8 inches of Double Layered Reinforced Concrete. 
My theoretical analysis in aerodynamics predicts in-flight stability with a minimum static 
margin of approximately two percent. The analysis and modeling of the electromagnetic 
launch resulted in an adequate design. For this analysis, I used three Comsol 
Multiphysics™ finite-element computer models. The modeling results validated 
fundamental railgun equations. The final projectile design concluded with a 3 m barrel 
and is characterized by the following parameters: conducting rails with an inductance 
gradient ~ 0.38 wtH/m; an average temperature rise in the rails of 20 °C per shot; an 


effective current of less than 2 MA; and a projectile launch velocity of 1100 m/s. 
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I. INTRODUCTION AND DESIGN OVERVIEW 


A. INTRODUCTION 

Several Department of Defense Development, Test and Evaluation (DT&E) 
activities along with other U.S. research organizations are exploring railgun technologies 
with the expectation of replacing current military munitions with hypervelocity 
projectiles. Of particular interest to the Army are currently funded programs to develop 
railgun applications for its Future Combat System (FCS) [1]. A key unresolved challenge 
in this effort is the design of an appropriate railgun projectile. The current primary 
concept for a potential FCS Army railgun projectile integrates a hypervelocity long-rod 
penetrator with an aluminum fall-away electrical conducting armature [2]. This projectile 
concept endeavors to meet the FCS requirement for defeating modern composite armor in 
conventional warfare. However, it has limited utility in urban military operations where 
the forward placement of dismounted infantry and minimization of collateral damage are 
important operational requirements. Given the current integration of units into combined 
arms teams and the likelihood of future unconventional urban conflicts, the optimal 
design of a railgun projectile must incorporate a multipurpose penetrator that minimizes 
fratricide and collateral damage. 

One objective of the FCS development program is to replace the capabilities of 
the M1 Abrams Main Battle Tank (MBT) and the M3 Bradley Infantry Fighting Vehicle 
(IFV) with one multipurpose armored vehicle [3]. The Army has already organized and 
fielded an Intermediate Brigade Combat Team (IBCT) that uses Light Armor Vehicles 
(LAV) in an attempt to replace the capabilities of the MBT and IFV. The IBCT still has 
two different vehicles to support combined arms operations, a mobile gun system and an 
infantry carrier vehicle. Major benefits of the LAV are that it is easily deployable via Air 
Force airlift assets and is highly mobile in rough terrain. The shortfalls of the LAV are a 
reduction in lethality and survivability [4]. The inherent deficiencies with the LAV have 
caused the Army to retain the MBT and IFV until it is able to field the FCS as more 
suitable replacement. 

The Army intends to develop a railgun for the FCS in order to match or exceed 


the lethality of the MBT and IFV, while trying to equal the deployability and mobility of 
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the LAV. The challenges for a railgun design include appropriate barrel dimensions, 
power supplies, and projectile launch configurations. Longer railgun barrels require less 
electrical current to achieve a hypervelocity projectile launch. The barrel length 
restriction for the FCS is approximately 3 meters. This estimate is based upon a barrel 
that is one-third the maximum vehicle length for the FCS which is eight meters [5]. A 
three meter railgun barrel is shorter than full scale experimental laboratory barrels (ten 
meters) [6]. The optimal geometries for the conducting rails are rectangular 
configurations. These geometries allows for the highest inductance per unit length, or 
inductance gradient, which directly affects the electromagnetic force on a projectile [7]. 
The primary candidate for the railgun power supply is a pair of inertial energy storage 
rotating machines, also known as Compensated Pulsed Alternators (CPA). The Army 
version of this power supply is projected to have a maximum energy storage capacity of 
20 MJ and a peak output current of 1.08 MA [8]. The Armor Piercing Fin Stabilized 
Sabot Discarding (APFSD) projectile is a leading candidate for the primary ballistics 
package because it has proven lethality against modern armored vehicles [2]. APFSD 
designers can easily configure the sabot petals to fit the geometry of the barrel [9]. The 
projectile launch package must include a device to conduct current and maintain 
electrical contact with the rails. The portion ofa railgun projectile’s launch package that 
conducts current induces the necessary electromagnetic force required for acceleration 
and final muzzle velocity. Army designers have decided that the best device for this 
purpose is a base-push or mid-ride fall-away aluminum armature. The proposed 
configurations of these armatures have shapes that allow the launch induced 
electromagnetic force to push the sides of this malleable metal against the rails [10]. 
Given the aforementioned FCS railgun constraints and design solutions, the Army 
believes that it can achieve launch velocities of 2.5 km/s using an APFSD projectile and 
deliver lethal impact energy to conventional armor targets. 

The Army designed the APFSD round for Kinetic Energy (KE) impacts against 
heavy to light skin enemy armored vehicles. The primary ballistic packages of APFSD 
projectiles are tungsten or depleted-uranium penetrators that have length to diameter 
ratios (L/D) of anywhere between 5 and 30, with respective diameters between 15 mm 


and 7 mm [Il]. These penetrators are highly effective at defeating modern enemy 
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armored vehicles. The high L/D characteristic that makes them devastating against armor 
actually work against them when penetrating urban targets. Bradley IFV gunners who 
have used conventional APFSD rounds against urban or concrete structures report that 
the penetrator passes through the target leaving a hole on the order of its body diameter. 
These gunners prefer to use the larger diameter aluminum-training round against concrete 
targets because it is significantly more effective [ll]. From a ballistic effects 
perspective, as seen by IFV gunners, the APFSD is not a multipurpose projectile. 

When considering engagements with the proposed railgun launch package, the 
fall-away armature and sabot petals induce restrictions on its safe battlefield employment. 
The U.S. military successfully trains with APFSD rounds and has used them in combat. 
The Department of the Army has mandated the following restrictions on the use of the 
APFSD with dismounted infantry: firing over the heads of unprotected personnel is 
prohibited; no personnel can be within 60 degrees of the gun target line for the first 
100 m of the projectile’s trajectory; no personnel can be within 7 degrees of the gun 
target line after the 100 m primary limit up to 400 m [12]. These restrictions are for the 
Bradley IFV using a conventionally fired APFSD, and they are even greater for the M1 
Abrams. The predictive error impact area, Surface Danger Zone (SDZ), for the Bradley 
M919 APFSD extends from the IFV out to 18 km [12]. The proposed railgun launch 
package will force the Army to impose more significant restrictions on the placement of 
dismounted infantry and SDZs because the aluminum fall-away armature and ballistic 
penetrator will have lethal launch velocities at 2.5 km/s. The Army cannot ignore these 
restrictions, and therefore, ground commanders will most likely relegate the use of 
railguns only to specific engagements such as those found in conventional armor warfare. 

If the urban and unconventional style of warfare in Operations Enduring Freedom 
and Iraqi Freedom are any indication of the future, then the proposed design of the Army 
railgun projectile will have limited utility. Armor and Infantry units operate in a 
combined arms team where dismounted infantry are often forward of supporting heavy 
armor assets. Within urban areas this is especially true. Armored vehicles that travel 
through urban terrain, without the support of forward or flanking infantry, are susceptible 
to side attacks by an unconventional enemy force. One must expect that dismounted 


infantry will continue to operate forward of FCS vehicles in future combat situations. 
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Another important issue in urban warfare is irregular friendly unit boundaries. Senior 
ground commanders often break up the control of towns or cities into multiple unit 
sectors. Subordinate commanders operating in these sectors normally find that they have 
weapons orientations toward other friendly units, and they must be mindful of the 
maximum effective range of their most lethal weapons. Given that company level sectors 
in urban operations are anywhere between 5 to 10 km in width or depth, the approximate 
18 km maximum effective range of a conventional APFSD will most certainly cross 
friendly unit boundaries. | Another significant consideration in urban operations is the 
minimization of collateral damage. Ground commanders cannot control or predict the 
amount of collateral damage that APFSD rounds will produce in urban terrain and are 
unwilling to use them unless there is a military necessity. Ifthe future FCS vehicle has a 
railgun with the proposed Army APFSD projectile, then firing safety issues will most 
likely constrain its employment [13]. 

In order for railgun technology to be more relevant to Armor and Infantry 
combined arms teams, the design of the projectile must incorporate the following 
concepts: 

e The projectile must be effective against various infantry targets, while 
minimizing collateral damage in urban warfare environments. 

e The projectile must be aerodynamically stable upon launch to provide 
precise trajectories for surgical strikes. 

e The projectile must be self-conducting during electromagnetic launch such 
that no parasitic mass falls away in flight, which minimizes SDZ and 


prevents possible fratricide. 


B. DESIGN OVERVIEW 

The primary objective of this thesis research was to create a credible and relevant 
alternative to the Army’s APFSD railgun projectile, i.e. a projectile that better satisfies 
the operational requirements of Infantry combined arms teams in future combat urban 
warfare scenarios. In order to prove the merits of a credible design alternative, this 
research involved the analysis of ballistic penetration, aerodynamic stability, and 


electromagnetic launch. I limited the scope of this research to theoretical design and 
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modeling. Using a recursive design loop involving ballistics, aerodynamics, and 
electromagnetic launch, I developed a plausible design alternative, which I refer to in this 


thesis as the Infantry Railgun Projectile. 


A critical step in any technical development process is to first create a set of 
reasonable constraints. Designers can then use these constraints to produce a baseline 
product for further testing and evaluation. After analyzing the limiting factors and 
operational requirements of an Army FCS railgun projectile, I created a specific set of 
constraints for my design (See Appendix A). Using these constraints, I developed a 
baseline 30mm primary diameter, 1.5 kg, design projectile that is shown in Figure |. I 
also created a scaled 40mm, 3.5 kg, projectile that I use to expose the advantages or 
disadvantages of size and weight. Detailed material properties, geometrics specifications, 
volume, mass and overall center of mass of both design projectiles are listed in the 


appendix (See Appendix B) [14], [15]. 


.«—Aluminum Armature 


” 


Flight Configuration 


PBlinsulator (Split Sections) — Tjtaniufn Fins 


a ee 
wm ore. 
‘ aa 


Figure 1. Infantry Railgun Projectile 


As shown in Figure 1, the Infantry Railgun Projectile is comprised of four major 
components. The tungsten penetrator serves as the primary ballistics package in the 
projectile. Its configuration allows for the successful penetration of armor targets and the 
perforation of concrete targets. The tungsten penetrator also acts as a secondary 
conductive armature to distribute the flow of current through the projectile. The 
aluminum armature is the principal electrical conductor in the projectile. It creates the 
first cross contact surface between the two conducting rails inside an electromagnetic 
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gun. The inverted “U” shape, cut into the base of the aluminum armature, allows the 
launch induced electromagnetic force to push its back edges against the conducting rails. 
This prevents “Transition” which is a technical term that railgun designers use to describe 
the loss of solid electrical contact between the conducting rails [7]. Since the electrical 
conductivity of tungsten is similar to that of aluminum, current will distribute itself 
through the tungsten penetrator and aluminum armature. The Polybenzimidazole (PBI) 
electrical insulator behind the aluminum armature prevents current from finding a cross 
conduction path between the rails. The insulator’s secondary purpose is to support and 
protect the deployable stability fins that are folded inside the projectile’s body during 
launch. PBI is a synthetic polymer that is temperature resistant up to 500 °C and is 
mechanically strong (42 MPa compressive yield and 160 MPa ultimate tensile strength) 
[15]. Another important property of this polymer is that manufacturers can easily cast 
this material into cylindrical shapes [16]. The last major component of the projectile is 
the titanium stability fins. Their function is to provide in-flight stability. The forward set 
of fins is for yaw stability, and the rear set is for pitch stability. I do not discuss the exact 
assembly of these components because that would involve specific engineering practices 
that do not fall into the scope of this research. Considering the sophistication of several 
Army weapon systems, I assume that the assembly of this projectile is both realistic and 
achievable. Overall, these four components act in concert with each other to allow for a 
successful electromagnetic launch, stable trajectory, and penetration of realistic targets. 
The configuration of the design projectile is suitable for electromagnetic launch 
from a rectangular railgun bore. The choice of materials and the geometric dimensions of 
the conducting rails are crucial factors in the analysis of electromagnetic launch. 
Multiple research organization have conceptualized, developed, and experimentally 
tested railgun barrels that use copper rectangular conducting rails. Given a generic 
version of this railgun barrel from the Institute for Advanced Technology (IAT), I have 
paired it with the design projectile to create a conceptual barrel that includes insulation 


and confinement [8]. The configuration of this barrel is shown in Figure 2. 


Infantry Railgun Barrel 
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Figure 2. Generic Railgun Barrel 


In summary, I developed both a baseline and a scaled projectile design that 
included constraints imposed by operational concerns, rectangular conducting rails and 
in-flight stability. I analyzed ballistic effects with semi-empirical equations. The ballistic 
analysis produced minimum penetration velocities that I feed into my aerodynamic 
analysis. Using semi-empirical aerodynamic equations, I made estimates on aerodynamic 
stability and found that the projectiles were at least neutrally stable at a flight velocity of 
1100 m/s with a minimum static margin of two percent. I then analyzed aerodynamic 
drag with further semi-empirical equations, and developed an equation for both 
projectiles that estimated their velocities with respect to flight distance. These equations 
allowed for calculations of final impact velocities for a 2 km and 500 m target. Using the 
impact velocities and a modified form of both design projectiles, I successfully modeled 
their impacts against the 2 km and 500 m targets with AUTODYN™, a finite-element 
hydro-code. Given the constraints of a 3 meter barrel length, an average temperature rise 
in the rails of 20 °C, and an 1100 m/s launch velocity, I successfully optimized the 


electromagnetic launch to produce final rail dimensions, an effective current, and the 
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inductance gradient. I used another finite-element code, Comsol Multiphysics™, to 
model static currents, magnetic fields, and the resulting electromagnetic force. I 
compared the modeled inductance gradients to their theoretic values, and found 
agreement to within 10 percent. Using this finite-element code one more time, I created a 
transient conductive heat model to analyze the temperature rise in the aluminum armature 
over the 5 ms launch time. I found that the aluminum armatures for both projectiles 
were below their melting points. In conclusion, the final designs of both Infantry Railgun 
Projectiles have satisfied the objectives of this thesis research. Although further work on 
these projectiles is necessary for their practical development, they serve as feasible 


alternatives to the APFSD railgun projectile. 


Il. BALLISTICS 


A. BALLISTICS OVERVIEW 

During the initial design phases of this research, I first tested all geometric and 
material projectile configurations against a primary ballistic objective. My ballistic 
analysis was always the first step in a recursive design loop that also analyzed 
aerodynamic stability and electromagnetic launch. I placed great emphasis upon 
choosing appropriate targets and evaluating theoretical penetration effects. The results of 
this ballistic analysis show that the design projectiles in this thesis have characteristics 


necessary to defeat modern battlefield targets. 


The primary ballistic objective for this research is the penetration of 100 mm of 
Rolled Homogenous Armor (RHA) with a zero degree obliquity angle (normal incidence) 
at a range of 2 km. The parameters of RHA thickness, obliquity angle and range are 
mutually important variables that represent an equivalent medium-armor target. All 
modern Infantry Fighting Vehicles (IFV) have kinetic energy projectiles. A comparison 
of international IFV armor penetrators shows diverse penetration capabilities [17]. 
Individual countries place different design importance in areas such as vehicle weight, 
rates of fires, ammunition carrying capacity, and survivability [18]. After investigating 
various IFV penetrators and their respective penetration capabilities, I found that the 
primary ballistic objective listed above is realistic. Furthermore, an IFV projectile that 


can penetrate this target is at the forefront of current technology [17]. 


I selected a standard 200 mm thick Double Layered Reinforced Concrete (DLRC) 
target as a secondary ballistic objective in order to evaluate the multipurpose 
characteristics of the railgun projectile. All modern IFV armor penetrators can easily 
perforate this target [11]. Army weapons researchers normally cite the perforation hole- 
diameter as an important ballistic effect in concrete targets [19]. Most IFV armor- 
piercing projectiles are long-rod penetrators. They are less than 20 mm in diameter, and 
therefore, leave small holes in DLRC. A modern IFV can fire multiple bursts at the same 
point of aim to produce a greater penetration diameter [11]. This tactic is effective 


against a thin non-reinforced concrete structure. When IFV gunners attempt this action 


against DLRC with armor piercing projectiles, they normally create a “jailed-window” 
effect in the concrete, shaped by the steel reinforcements [11]. In certain circumstances, 
ground combat commanders want to use an IFV to create breech-holes through buildings 
and reinforced bunkers [11]. For this reason, researchers want to develop IFV 
penetrators that will create the largest possible hole in concrete. Several experts agree 
that a 200 mm DLRC target is a realistic choice when the focus of research is to analyze 


penetration and perforation effects against urban concrete structures [19]. 


Throughout the 20" century, ballistic experts used both empirical and semi- 
empirical formulas to evaluate the penetration capability of solid cylindrical projectiles 
against armor targets [20]. Before the use of computers, some modern experts focused on 
the limit velocity to gage penetration effects. Researchers developed and refined a semi- 
empirical formula, the velocity ballistic limit, to estimate minimum penetration velocities 
in thick armor targets. Ballistic experts restrict the application of this formula to specific 
projectile geometries, armor thicknesses, and impact velocities. The parameters of the 
projectiles designed in this thesis satisfy the constraints for the formula and match the 
range of experimental data that researchers used for its validation. I compared the 
accuracy of the velocity ballistic limit against one set of equivalent experimental data 
[21]. The velocity ballistic limit agreed with the results of the external experimental data 
to within two percent. In this report, I initially used the velocity ballistic limit formula as 
the primary analytical method for analyzing RHA penetration effectiveness. The major 
utility of this formula is that it allows one to calculate the limit velocity or the minimum 
impact velocity required to create penetration in a particular armor target. I needed this 
estimate in order to conduct my analysis for aerodynamic stability and electromagnetic 


launch. 


Ballistic experts turned their interest toward reinforced-concrete impact during the 
mid 1940s [22]. By 1946, the National Defense Research Committee (NDRC) developed 
a penetration, scabbing, and perforation thickness empirical formula. In 1966, 
researchers adjusted these formulas to include the compressive strength of concrete. 
They refer to the new formulas as the modified-NDRC [22]. As recent as 2002, a study 
in the International Journal of Impact Engineering validated the use of these formulas for 


velocities up to | km/s [23]. Ballistic experts restrict the application of these formulas to 
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specific projectile geometries, concrete thickness to projectile diameter ratios, and impact 
velocities. The parameters in this thesis satisfy the constraints for the formulas and 
match the range of experimental data that researchers used for its validation. I used the 
modified-NDRC formulas as the primary analytical methods for estimating penetration, 


scabbing, and perforation effects in the 200 mm DLRC target. 


In the absence of experimental testing, computer finite element modeling serves 
as the next best refinement to analytical estimates. For the purpose of validating ballistic 
estimates, I used the AUTODYN finite-element computer model to simulate impacts 
against RHA and DLRC. The modeling domain was two-dimensional (2D) axial- 
symmetric. I accounted for my design projectile’s non-uniform diameter by adjusting 
this parameter to reflect an equivalent value. I reorganized the orthogonal steel- 
reinforcement structure in DLRC into circular rings in order to satisfy the axial- 
symmetric constraint in AUTODYN. My comparison of the RHA computer model to the 
velocity ballistic limit resulted in a greater penetration capability for the computer model. 
The AUTODYN stimulation for impacts into DLRC validated the modified-NDRC 
estimates for penetration, scabbing, and perforation. AUTODYN ballistic simulations 
are extremely useful in representing near realistic projectile impact scenarios. However, 
they cannot stand alone without agreement between theoretical estimates or experimental 


verification. 


In order to better understand the discrepancy between the RHA computer model 
and the velocity ballistic limit formula, I concluded my analysis with a refined theoretical 
analysis. Using a variant analytical method of the Tate Model, which researchers 
commonly refer to as critical analysis, I found agreement with the RHA computer model. 
In conclusion, I am confident in the penetration effectiveness of the design projectiles 
against my primary ballistic objective, RHA. My analysis of the perforation effects in 
DLRC indicates that the design projectiles will create holes in this target as large as their 
primary diameters. 

B. MODIFIED PROJECTILES 

I created two modifications to the original configuration of the design projectile. 

These modifications were necessary to conduct theoretical estimates and computer model 


simulations. The first modification, MOD-1, pertains to calculations with the velocity 
ll 


ballistic limit and the modified-NDRC formulas. The second modification, MOD-2, 
pertains to the AUTODYN finite element analysis. Both modifications are conservative 
adjustments to the original design and involve reasonable assumptions in the choice of 


each characteristic parameter for the two modifications. 


When researchers use the velocity ballistic limit or the modified-NDRC formulas 
to evaluate ballistic effects, they can only conduct analysis on projectiles that have 
circular-cylindrical geometries with conical nose cones. Given this shape, the formulas at 
most require the characteristic diameter, length, mass, and nose fineness. If we consider 
the geometry and material composition of the design projectile in this thesis, we can 
assume that the tungsten portion of the projectile serves as the primary ballistic 
penetrator. Therefore, it is the only shape that one needs to consider in a ballistic 
analysis. The sides of the tungsten penetrator are shaved, and are the only geometric 
configurations that prevent it from being a perfect circular cylinder. In order to analyze 
the penetrator as a perfect cylinder, I modified the geometry to reflect a characteristic 


diameter, length and mass. I refer to this modified projectile as MOD-1. 


It is important to define the characteristic diameter, length, and mass of MOD-1 in 
order to prevent any confusion in this ballistic analysis. The characteristic diameter is the 
equivalent diameter, d, of the tungsten penetrator’s major body cylinder. One can 


calculate the equivalent diameter by assuming that the cross-sectional surface area, S,,;, 


is equal to the area of a perfect circle, S,,., =(ad°/4). The result of this calculation 


produces equivalent diameters of 29 mm and 38.8 mm for the respective 30mm and 
40mm design projectiles (See Chapter III, Section C.1.b.). I assumed that the 
characteristic length of MOD-1 is the distance from the tip of the nose to the end of the 
major body cylinder. This is a reasonable assumption since both the nose cone and major 
body cylinder are the only geometries that will force a perforation hole through an armor 
or concrete target. The characteristic mass of MOD-1 is the mass of the complete 
projectile, which includes materials that have no ballistic relevance. One must consider 


the total mass because both the velocity ballistic limit and the modified-NDRC formulas 


use relationships between impact kinetic energy, MV* / 2, and displaced target volume. 
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One can further examine the general shape of MOD-1 and the exact specifications for it 


with regard to both the 30mm and 40mm design projectiles in Figure 3. 
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Figure 3. Projectile Modification for Theoretical Analysis (MOD-1) 


The AUTODYN modeling domains that I used in my computer simulations were 
two-dimension (2D) axial-symmetric. This modeling domain allowed the input of 
irregular 2D geometries. Using the same equivalent diameter from MOD-1, I developed 
a second modification for the AUTODYN analysis. With the ability to include more of 
the original projectile, I included the tungsten penetrator’s rod cylinder up to the leading 
edge of the first tail fin slot, because the fin slot cuts directly through the center of the rod 
cylinder. The effective mass of MOD-2 is based upon the density of tungsten and the 
physical geometry of the modified projectile. MOD-2 has a mass that is 82 percent of the 
original design projectile (See Appendix B). Although the mass of MOD-2 is a 
significantly conservative value, the geometric constraints that I applied were necessary 
given the near realistic environment of the AUTODYN analysis. The general shape and 


parameter specification of MOD-2 are given in Figure 4. 
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Properties of Tungsten 
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Figure 4. Projectile Modification for Computer Modeling (MOD-2) 


C. TARGETS 

i 100 mm of Rolled Homogenous Armor 

Both U.S. and international military ballistic laboratories use RHA as the standard 
by which to compare armor protection or ballistic penetration. In a time when new 
composite armor is multi-layered and complex, weapons researchers have found it useful 
to rate this armor in an equivalent RHA thickness. When rating a target, one must state 
the angle of obliquity, range to target, and RHA thickness. At a fixed target range, if we 
rate a piece of armor at an equivalent RHA thickness of 50 mm with a 60 degree angle of 
obliquity, we can approximately gage that same armor with an 85 mm thickness at 
normal incidence. One can understand the validity of this statement after examining the 
effects of the obliquity angle in the velocity ballistic limit formula, which I develop later 


in this chapter. 


In an effort to select the appropriate RHA thickness for an IFV armor penetrator, I 
investigated the ballistic effectiveness of several medium-caliber Armored Fighting 
Vehicle (AFV) munitions. My information source, “Jane’s International Defense 
Review,” listed at least two modern 25mm-caliber projectiles with an RHA penetration 
capability of 31 mm. For both projectiles, the author listed an impact obliquity angle of 
60 degrees and target range of 2 km. The author identified one modern 30mm-caliber 


projectile with an RHA penetration capability of 47 mm at an obliquity angle of 60 
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degrees and target range of | km. For that same caliber, the author also referenced a state 
of the art projectile with a 120 mm RHA penetration capability at normal incidence, 
given a target range of | km [17]. Using this information, I selected an RHA thickness of 
100 mm, a zero degree or normal incidence obliquity angle, and a range to target of 2 km. 
I chose these target parameters in order to make the design of the Infantry Railgun 
Projectile competitive against other medium caliber penetrators. Figure 5 shows the 
target and its material properties. I made the diameter of this target sufficiently large 
such that the cross-sectional area of the impacting design projectile was less than 1% of 
the target area. This minimizes any radial expansion of the target during impact and 


properly reflects the surface area of a large armored vehicle. 


Rolled Homogeneous Armor (RHA) 
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Figure 5. RHA Target Diagram 


23 200 mm of Double Layered Reinforced Concrete 

An Army Research Laboratory (ARL) report on hypervelocity penetration into 
concrete targets motivated my selection of a DLRC target. The ARL report stated that its 
DLRC target represents a typical concrete target in a modern urban battlefield [19]. The 


specifications of the Full Scale DLRC target are shown in Figure 6. 
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Figure 6. Full Scale DRLC Target Diagram 


The orthogonal reinforcement geometry in the DLRC target is not suited for the 
AUTODYN 2D axial-symmetric modeling domain. After consultation with my primary 
ballistics advisor, I modified the orthogonal structure into circular reinforcements [24]. I 
also decided to change the reinforcement material from Grade 40 to 4340 Steel. The 
motivation for this change came from the absence of Grade 40 steel as a selectable 
material in the AUTODYN Material’s Library [25]. I increase the compressive strength 
of the DLRC to 34 MPa because that value better represents a concrete barrier that is 
designed to resists projectile penetration. As with the RHA target, I made the diameter 
of the modified DLRC target significantly larger than the diameter of the impacting 
design projectile. This minimizes any radial expansion of the simulated target and 


properly reflects a larger wall surface. The modified target 1s shown in Figure 7. 
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Figure 7. Modified DLRC Target Diagram 

I selected a maximum range to target of 500 m in order to represent a realistic 
engagement situation within an urban environment. If we compare the impact velocities 
that I used in the RHA analysis against those used for DLRC then we will find larger 
velocities for the concrete impacts. This difference comes from my estimate of impact 
velocities. I made these estimates by using an aerodynamic drag analysis, which | 
discuss later in this report (See Chapter III, Section F). 
D. THEORETICAL PREDICTIONS FOR IMPACTS INTO RHA 


I used the velocity ballistic limit model to estimate minimum penetration 
velocities for both design projectiles into 100 mm of RHA. My estimates for the 
maximum allowable impact velocities are based upon constraints from my aerodynamic 
drag analysis (See Chapter III, Section F, and Table 9). The existence of known impact 
velocities enables one to compare these values against the ballistic limit and to 
simultaneously evaluate residual velocities after impact. The velocity ballistic limit 
model predicts that only the 40mm design projectile can penetrate the 100mm RHA 
target. The success of at least one of the proposed projectiles demonstrates a need for 


validation using comparable theoretical techniques or computer model simulations. 


A scientist by the name of de Marre proposed the first version of the ballistic 
limit formula in 1886. He assumed that there was a physical relationship between the 


limit energy of an impacting projectile and the amount of target volume that it displaces. 
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The limit energy is the minimum kinetic energy of a penetrating projectile [20]. 
Although the formula that de Marre proposed was dimensionally awkward, its underlying 
assumption is that twice the limit energy is proportional to the perforation volume in the 
target. We can better examine this relationship in the following expression: MV* « 7D” ; 
where is the mass of the projectile, V is the limit velocity, 7 is the target thickness, and 
D is the projectile’s characteristic diameter [20]. Throughout the early 20" century, 
several theorists and ballistics experts refined the original formula into the modern 
velocity ballistic limit model. They also used a semi-empirical momentum relationship 
between the penetrator and the ejected target material to devise a residual velocity 
relationship [20]. The residual velocity relationship makes estimates from the input 
parameters and the results of the velocity ballistic limit. Figure 8 shows a summary of 


the entire velocity ballistic limit model. 
Velocity Ballistic Limit Residual Velocity 
0:[7, <¥] 
ro ip. 
a? -V?)"? [V, > ¥ ] 


Description and Units 


Vi = limit velocity [m/s] 
V+ = residual velocity [m/s ] 
« = (4000)"2 for RHA targets 
£= length of penetrator [ cm ] 
D = diameter of penetrator [ cm ] 


M'= PrAD *z/A4 ¢ = 0.3; for long rod penetrators 
M k= 0.75; for long rod penetrafors 
pp I= mass of penetrator [ g ] 
M+M'3 T = target thickness [ cm ] 


#= obliquity angle [ deg ] 
Mf’ = projected impact mass of target, 
prior to impact [ g ] 
pT = density of target, 
7.85 for RHA [ g / ce] 


p=24+2/3 


Figure 8. Velocity Ballistic Limit Model 
There are several constraints that one must consider in order to properly use the 
velocity ballistic limit. The parameters c and k, listed in Figure 8, are only well defined 


for specific long-rod penetrators and target thicknesses. A researcher by the name of 
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Lambert empirically derived these parameters from an experimental data base of 200 
impacts involving penetrators with the following characteristics: mass between 0.5 and 
3630 grams; Length to Diameter (L/D) ratio between 4 and 30; impact obliquity angle 
between 0 and 60 degrees, and material density between 7.8 and 19 g/cc. The target 
thickness in Lambert’s limit-velocity database was between 0.6 and 15 cm [20]. Ballistic 
experts who have analyzed Lambert’s parameters against his database believe that the 
ratio of target thickness, 7, to penetrator diameter, D, must satisfy the following 
relationship: 7/D > 1.5. They have further added the following general constraints to 
the geometry and configuration of the penetrator. It should have a right-circular 
cylindrical body with a conical, hemi-spherical, or flat nose. There should be obvious 
values for the diameter, length and mass. If geometries are approximately cylindrical 
then one can use effective values to estimate the characteristic diameter, length and mass 
[20]. The 30mm and 40mm design projectiles and the 100mm RHA target fit all of the 
constraints listed above, and one should have confidence in the reliability of limit- 


velocity estimates for this design. 


In order to gain a higher level of confidence in the velocity ballistic limit, I tested 
the formula against one set of external data that also fit within the aforementioned 
restrictions for penetrator and target dimensions. This experimental data comes from the 
1971 Graberek database [21]. I searched for an experimental data point that gave the best 
match to the type of RHA impact in my design. The closest data point was a 20 gram 
tungsten penetrator, with a 4 mm diameter and an L/D ratio of 20, impacting a 50mm 
RHA target. The experimental impact was at normal incidence, and the reported limit 
velocity was 1260 m/s. Given these specifications and results, I solved explicitly for the 
parameter a, which for RHA should be (4000)°. My result for a is(4057)°. This gives 
a percent error of 1.4 %. There were no other similar data points in Garabrek’s 
database. I assume that this agreement with external data further validates the reliability 


of the velocity ballistic limit formula. 


I applied the velocity ballistic limit model to an impact involving the modified 
penetrator, MOD-1, and the 100mm RHA target. Given the input parameters for MOD-1 
and the RHA target, I calculated limit velocities of 1028 m/s and 844 m/s for the 
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respective 30mm and 40mm design projectiles. My aerodynamic analysis imposed a 
launch velocity constraint of 1100 m/s (See Chapter III, Section D). The impact 
velocities for the 30mm and 40mm projectiles based upon in-flight velocity reduction due 
to drag is discussed in the aerodynamics section (See Chapter II, Section F, and Table 9). 
An examination of these impact velocities shows that the 30mm projectile’s limit velocity 
is higher than the aerodynamically calculated impact velocity; whereas, the 40mm 
projectile’s limit velocity is lower. Therefore, the velocity ballistic limit model predicts 
that only the 40mm projectile will penetrate the RHA target. Using the impact velocity 
and the remaining portion of the model, I calculate a 413 m/s residual velocity after 
penetration for the 40mm projectile. The results from the velocity ballistic limit model 
require further validation through a comparable theoretical method or computer model 


simulation. The input parameters for the model and results are given in Table 1. 


Projectile | Projectile 
| M_|_g_ |MassofMOD-1_————sd|sCs1497_—«| 3549 
| @ | deg | ObliquityImpactAngle = | 0 =| 0 Tf 


ore VL & Vr Parameter 3.45 2.57 

= VL Parameter_ 2.48 1.65 

= ee ee Impact Veloci 828 890 
|g | MassofEjected Target | 
ee 


Beet Vr Parameter 3.15 2.86 


G m/s | Calculated Residual Velocity 0 413 
Table 1. | Penetration Estimates for Projectile Impact into RHA 


E. THEORETICAL PREDICTIONS FOR IMPACTS INTO DLRC 

The modified-NDRC’s concrete penetration, scabbing, and perforation formulas 
are the primary analytical methods that I use to estimate ballistic effects in the DLRC 
target. The results of the NDRC formulas predict penetration, scabbing and perforation 
for impacts involving both design projectiles and the 200mm DLRC target. 
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In the last half of the 20 century, before the use of finite element computer 
codes, ballistic researchers developed empirical formulas to estimate military ballistic 
effects against reinforced concrete. There are multiple empirical formulas that estimate 
penetration depths and perforation diameters from the normal impact of solid cylindrical 
projectiles against light to moderately reinforced concrete targets [22]. Researchers have 
conducted critical comparisons of the penetration formulas as recent as 2001. A research 
report in the International Journal of Impact Engineering states that the modified-NDRC 
penetration formula gives the best agreement with experimental data over a broad range 
of projectile geometries and impact velocities up to 1 km/s [23]. American military 
researchers originally developed the NDRC scabbing and perforation formulas from an 
experimental database that included high velocity impacts of steel cylindrical missiles, 
with diameters of between 39 mm and 155 mm into reinforced concrete with 7/D ratios 
between 3 and 18 [22]. Ballistics experts restrict the use of the NDRC formulas to the 
following impact parameters: impact velocities greater than 152 m/s, the impacting 
projectile is non-rotating, and impacts are at normal incidence. The proposed impacts of 
the 30mm and 40mm railgun projectiles into the 200mm DLRC target approximately fit 
the parameter range that the NDRC used to develop its formulas and also meet the impact 
parameter restrictions just mentioned. Therefore, one should have confidence in any 
predictive estimates from the application of the NDRC formulas in this design. The 


modified-NDRC formulas for penetration, scabbing and perforation are show in Figure 9. 


Modified-NDRC Penetration Description and Units 
130NW V x = penetration depth [inch ] 
x=] & = projectile diameter [inch ] 
ale 1000d N = nose shape factor; 
= 1.4 for very sharp noses 
Modified-NDRC Scabbing W = projectile weight [ Ibs ] 
x ¥ = impact velocity [ ft/s ] 
{0653 5133} Je= concrete compressive 
strength [ psi] 
Modified-NDRC Perforation és = concrete thickness to resist 


t, =1.32d+1.24x {135s4.<135h scabbing [ inch ] 
P a tp = concrete thickness to resist 


perforation [ inch ] 


Figure 9. | Modified-NDRC Formulas (British Units) 
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American military researchers developed the NDRC formulas in order to estimate 
required concrete thickness to prevent penetration, scabbing and perforation from 
ballistic missiles. I use these formulas to show that the proposed 30mm and 40mm 
railgun projectiles can penetrate and produce scabbing and perforation on the 200mm 
DLRC target. If we apply the input parameters from MOD-1 projectile, for impact 
velocities that I calculated in my aerodynamic analysis, then we can calculate the 
required thickness of the DLRC target to prevent penetration, scabbing, and perforation. 
Table 2 shows the input parameters for the NDRC formulas and the resulting concrete 
thickness needed to prevent penetration, scabbing and perforation. All of the calculated 
concrete thicknesses are at least eight times the thickness of the 8 inch or 200 mm DLRC 
target. The required thickness to prevent scabbing and perforation are larger than the 
required thickness to prevent penetration. The expectation of needing larger concrete 
thickness to prevent scabbing or perforation, given that the projectile does not penetrate, 
comes from shock wave and shear concrete plug pushed forward by the impacting 
projectile. The combination of the shock wave and concrete plug can create scabbing and 
perforation even though a projectile does not completely penetrate. The results of the 
NDRC formulas predict that both design projectiles will penetrate and produce scabbing 


and perforation on the DLRC target. 


Variable 


Description 30 mm Projectile 40 mm Projectile 


Metric 


British 


Metric 


British 


Diameter of MOD-1 


29 mm 


1.14 in 


38.8 mm 


1.53 in 


Impact Velocity 
Weight of MOD-1 


1024 m/s 
1497 g 


3359 ft/s 
3.29 lbs 


1043 m/s 
3549 g 


3421 ft/s 
7.81 Ibs 


d 

Vy 
W 
N 


Nose Shape Factor 


1.14 


1.14 


1.14 


1.14 


Compressive 
Strength DLRC 


34 MPa 


4900 psi 


34 MPa 


4900 psi 


DLRC Thickness to 
Prevent Penetration 


1742 mm 


68.6 in 


2514 mm 


99.0 in 


DLRC Thickness to 
Prevent Scabbing 


2431 mm 


95.7 in 


3502 mm 


137.9 in 


DLRC Thickness to 
Prevent Perforation 
Table 2. 


2199 mm 


ae 


86.6 in 


3169 mm 


Input Parameters and NDRC Results 


124.8 in 


Although reinforcement does not explicitly appear as an input parameter in the 
NDRC formulas, these formulas inherently consider concrete reinforcement in their 
calculations. The NDRC developed its formulas from experimental data using non- 
deformable solid cylindrical projectiles against light to moderate reinforced concrete 
targets. The modified-NDRC goes one step closer to a semi-empirical formula by 
including the compressive strength of concrete, fc, which the original formula did not do 
until its modification in 1966. A researcher that is associated with the American Society 
of Civil Engineers (ASCE) evaluated the NDRC formulas for their lack of specificity 
with regard to the amount of concrete reinforcement. He concluded, from the analysis of 
experimental data, that the amount of reinforcement is a negligible concern unless the 
projectile diameter is equivalent to the reinforcement bar spacing [22]. Ballistic experts 
generally accept the fact that steel reinforcement within concrete serves to maintain 
structural integrity after impacts or explosions. Heavy reinforcement may restrict the 
perforation hole to the diameter of the impacting projectile. Reinforcement bars are not 
designed to prevent penetration [19]. Both design projectiles have significantly smaller 
diameters than the bar spacing of the DLRC target. From the conclusions of the report 
from the ASCE, one should be able to neglect the amount of reinforcement when 
considering NDRC estimates in this design. The fact that the 200mm DLRC target has 
heavy reinforcement suggests that the impacting 30mm and 40mm design projectiles will 
create holes in the concrete that are slightly larger than their primary diameters. 

F, BALLISTIC MODELING 

In order to refine or validate my theoretical estimates, I simulated the design 
projectiles impacting into the RHA and DLRC targets at their calculated impact 
velocities. I used the AUTODYN version 5.0 interactive non-linear dynamic analysis 
software for all simulations [25]. The AUTODYN program offers planar or axial 
symmetric coordinates in two or three dimensions. It also allows the user to choose 
between Eulerian or Lagrangian reference frames, each with their own advantages or 
disadvantages to a particular model. The equation of state that a user chooses to represent 
materials in the modeling domain is an important consideration. Given the sub- 
hydrodynamic impact scenarios of solid cylindrical tungsten projectiles approaching their 


solid targets at normal incidence, with expected deformation of both projectile and target 


ZS 


during impact, I choose to operate the AUTODYN simulation in a Lagrange, 2D axial 


symmetric coordinate system with linear equations of state for materials. 


The normal incidence direction of impact and the inherent axial symmetry of the 
projectile allows for 2D axial-symmetry. The only problem with this coordinate system, 
as it applies to the proposed impact scenarios, is for the original DLRC target. The 
modification of the DLRC target into an axial symmetry representation was a necessary 
but reasonable change. AUTODYN allows users to map the 2D axial symmetric solution 
onto a circular three dimensional (3D) representation of the impact scenario. Although 
this is a computational shortcut, it adequately represents the 3D interactions of materials 


in computer simulation models. 


The Lagrange reference frame follows the interaction of each mass particle in the 
modeling simulation. The Euler reference frame calculates the flow of mass through 
static cells. It is important to follow the interaction of individual mass particles for sub- 
hydrodynamic or less than 3 km/s impacts. In this impact regime, mass particles remain 
solid and respond to interaction based upon elastic or plastic response patterns. Solid 
materials respond in a plastic or deformable manner within the shock front created by a 
high velocity impact. The AUTODYN equations that calculate interaction, based upon 
the Lagrange reference frame, ensure the conservation of mass, momentum, and energy. 
AUTODYN’s ability to ensure these conditions are met, while allowing each particle to 


deform, creates the most realistic modeling simulation for this design. 


When a mass particle is allowed to deform, it is important to consider the zoning 
that one applies to mass particles. The amount of zoning that one chooses directly 
influences computational memory and the time steps taken by a computer simulation. In 
all simulations, the zoning was set to two cells per unit millimeter. Therefore, one square 
millimeter had four zones. This zoning was applied to the entire model except for the 
nose cone of the projectile. The base of the projectile was set to two zones per unit 
millimeter. However, the zoning increased towards the tip of the nose. The type of 
geometry that AUTODYN allowed for a triangular configuration set the structure of the 


zoning in the nose. Given the high level of deformation that a nose cone experiences 
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during impact, the AUTODYN zoning in the projectile’s nose was both reasonable and 


necessary. 


The equation of state that one chooses for materials is also a selectable modeling 
characteristic, and needs to properly represent the range of interactions within a particular 
model. The equations of states for all materials were linear Von Mises, except for the 
concrete and the steel reinforcement bars. The linear Von Mises model evaluates mass 
particles for any violation of their material-specific ultimate yield stress. If the model 
finds a violation of the Von Mises yield criterion, it linearly adjusts cells to satisfy 
appropriate boundary conditions [25]. The proper equation of state for the concrete is 
one that represents a porous material. This is a selectable equation of state in 
AUTODYN [25]. The 4340 steel that I used to represent the reinforcement material had 
a default linear Johnson-Cook model. This model allows a simulation to calculate the 
effects of compression with temperature. Since the reinforcement bars did not interact 
directly with the impacting projectiles, their equation of state was not an important design 
parameter. The default Johnson-Cook model was adequate. Overall, the materials in 
RHA and DLRC simulations start out as solid and remain solid. Their response to 
pressure and physical interaction is either elastic or plastic based upon their proximity to 


the shock front created by the impact. 


When using the AUTODYN modeling software, users most consider the effects 
of reference frames and material interactions. I ensured that the impact scenarios 
between the tungsten penetrator and both the RHA and the DLRC targets were setup to 
reflect the proper interaction of mass particles. A summary of the program set-up for 
each impact scenario is listed in the appendix (See Appendix D). 

1. AUTODYN Simulation of MOD-2 Impacting into RHA 

I executed individual AUTODYN modeling simulations to predict the ballistic 
effects of the 30mm and 40mm design projectiles impacting into the 100mm RHA target. 
AUTODYN’s 2D axial symmetrical modeling domain required the use of projectile 
modification MOD-2 for each projectile. I did not change the specification of the 100mm 
RHA target that I previously identified in the Targets section. The 30mm and 40mm 
MOD-2 projectiles impacted the RHA target at their respective impact velocities of 828 


m/s and 890 m/s. The results of the modeling simulations predict that both MOD-2 
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projectiles will perforate the RHA target. The residual velocities of the 30mm and 40mm 


projectiles are 758 m/s and 847 m/s, respectively. 


The AUTODYN results for the 30mm projectile differ significantly from the 
estimate of the velocity ballistic limit, which predicts no penetration or residual velocity. 
The result for the 40mm projectile agrees with penetration prediction of the velocity 
ballistic limit. The AUTODYN residual velocity prediction for the 40mm projectile is 
over two times larger than the 413 m/s estimate of the velocity ballistic limit. The 
disagreement between the AUTODYN models and the velocity ballistic limit estimates 
requires further theoretical or experimental analysis. The velocity plots and final 
graphical images of the two impacts are shown in Figures 10-13, under the subsections 


for each projectile. 


a. 30mm MOD-2 Impact into RHA 
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Figure 10. Simulation Velocity of 30mm Projectile During Impact with RHA Target 
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Figure 11. 


40mm MOD-2 Impact into RHA 
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Figure 13. Simulation Image of Projectile and RHA Target 4 ms After Impact 


2. AUTODYN Simulations of MOD-2 Impacting into DLRC 

Using the AUTODYN finite element software, I simulated the individual impacts 
of the 30mm and 40mm design projectiles against the 200mm DLRC target. The 2D 
axial symmetrical modeling domain within AUTODYN required the use of projectile 
modification MOD-2 for each projectile. The original DLRC target that I previously 
specified in the Targets section required a modification into a proportionate 2D axial 
symmetric geometry. The modified DLRC target is also discussed and graphically 
represented in the Targets section. The 30mm and 40mm MOD-2 projectiles impacted 
the DLRC target at their respective impact velocities of 1024 m/s and 1043 m/s. The 
results of the modeling simulations predict that both projectiles will penetrate and 
produce perforation and scabbing in the DLRC target. The residual velocities of the 
30mm and 40mm projectiles are 997 m/s and 1023 m/s, respectively. However, the 
residual velocities have no relevance toward estimates from the modified-NDRC 
formulas. The graphical AUTODYN results match the predictions from the modified- 
NDRC formulas. The simulations do not show the formation or backside ejection of a 
concrete plug. The ratio of the length of the nose to the base diameter, or nose fineness, 


of each projectile is 3.4. This is a high value for nose-fineness and causes a higher lateral 
28 


distribution of the normal shock front. If the noses of the projectiles were blunt then we 
would expect to see the formation of a concrete plug. The minimum perforation 
thickness in each simulation was approximately 2 mm larger than each projectile’s body 
diameter. The NDRC formulas do not allow for a predictive estimate of this value. 
Therefore, one can only conclude that the minimum perforation holes created by these 
projectiles in DLRC will be no smaller than their primary diameters. The velocity plots 
and final graphical images of the two impacts are shown in Figures 14-17, under the 
subsections for each projectile. 


a. 30mm MOD-2 Impact into DLRC 
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Figure 14. Simulation Velocity of 30mm Projectile During Impact with DLRC 
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Simulation Velocity of 40mm Projectile During Impact with DLRC 


Figure 16. 
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Figure 17. Simulation Images of 40mm Projectile and DLRC Target 


G. CRITICAL ANALYSIS 

The disagreement between the RHA impact estimates of the velocity ballistic 
limit and the results of the AUTODYN simulation required that I make further 
investigations into either theory or modeling. A variant of the Tate Model, called critical 
analysis, fit the sub-hydrodynamic impact scenario of a rod penetrator traveling at less 
than 3 km/s into a solid metal target [26]. The theoretical foundations of the Tate Model 
come from the pressure balance equation at the stagnation point between an impacting 
projectile and the target, as shown in Figure 18. Furthermore, this model assumes that 
the elastic response regions of both the projectile and target material provide residual 
strength to counter deformation. A physical description of this equation is shown in 


Figure 18. 
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YP = Materiel Slrenylt of Projectile 
RT = Material Strength of Target 

Ps = Stagnation Pressure in Projectile 
PT = Stagnation Pressure in Target 


VP = Velocity of Projectile 

VT = Cavity Propagation Speed 
pP = Density of Projectile 

pT = Density of Target 


Figure 18. Sub-hydrodynamic Penetration Model 


The Tate Model assumes that the velocity of the projectile does not significantly 
change during penetration. Specifically, it assumes that difference between the velocity 
of the projectile and the cavity propagation speed is a constant for a penetrator that has a 
relatively small longitudinal shock wave velocity. The critical analysis method diverges 
from the Tate Model in that it assumes the projectile’s velocity changes during impact. It 
focuses calculation around a critical velocity, VC. This is the point in the impact where 
the velocity of the projectile, VP, equals the cavity propagation speed, VT. One can 
apply this method to two different impact scenarios in which the material strength of the 
projectile is either greater than or less than the material strength of the target. Applying 
the critical analysis method to the impact of a tungsten projectile, with a material strength 
of YP = 1.51 GPa, against a RHA target, with a material strength of RT = 1.32 GPa, the 
appropriate case is where the strength of the projectile is greater than the strength of the 
target, YP > RT [26], [27]. 
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Given the situation in which the strength of the projectile is greater than the 
strength of the target, one can solve for the cavity propagation speed, VT. The analytical 
expression for the cavity propagation speed is complicated, and I have decided to 
reference it and further important relationships in the appendix (See Appendix D). As I 
stated earlier, the projectile’s critical velocity is approximately the cavity propagation 
speed. One can develop a relationship for the penetration depth per length, P(VC)/L, of 
the projectile for the portion of the impact where the projectile’s velocity is greater than 
the critical velocity. When the projectile is less than the critical velocity, we can assume 
that it continues to penetrate the target as a rigid projectile with depth Prigid/L. For rigid 
projectile penetration, the projectile continues to slow down until it either passes through 
the target or reaches a velocity of zero. One can calculate the total penetration depth by 


summing the both components of the penetration per unit length of the projectile, 


BPC) Pie 


7 7 [26]. Again, the relationships and equation manipulations for this 


analysis are complicated and left to the appendix (See Appendix D). 


The results of the critical analysis for both the 30mm and 40mm design projectiles 
impacting the RHA target indicate a successful penetration for each impact scenario. The 
total penetration depth of the 30mm projectile is 1.24 times the projectile length. The 
total length of the nose and major body cylinder of the 30mm projectile is 150 mm. If 
one includes the rod cylinder then the total length is 261 mm. It is important to analyze 
just the nose and major body cylinder because these aspects of the projectile have the 
largest characteristic diameter. If they do not have enough length to penetrate the RHA 
target then it is likely that they will plug the impact path and not allow the smaller 
diameter rod cylinder to continue penetrating. Considering just the nose and major body 
cylinder of the projectile, the total penetration depth is estimated to be 186 mm, which is 
greater than the 100 mm thickness of the RHA target. The calculated penetration depth 
of the 40mm projectile is 1.3 times its length. Again, considering the total length of the 
nose and major body cylinder which is 200 mm, then one should calculate a total 
penetration depth of 260 mm. This is significantly larger than the thickness of the RHA 


target. 
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Although the critical analysis does not yield a comparative limit velocity, it does 
indicate that the projectiles have enough length to successfully penetrate the 100mm 
RHA target. An implied conclusion from the critical analysis is that both projectiles will 
have a significant residual velocity. One can only approximate the residual velocity by 
assuming that it is close to the cavity propagation speed. This value is 523 m/s for the 
30mm projectile and 562 m/s for the 40mm projectile. Assuming that the cavity 
propagation speeds represent the residual velocities from the critical analysis, we can 
compare them to the AUTODYN predictions of 758 m/s and 847 m/s for the respective 
30mm and 40mm design projectiles. The percent difference for the 30mm projectile is 
31% and for the 40mm projectile it is 33.5%. Overall, the critical analysis gives better 
agreement with the AUTODYN simulation for RHA impact, than does the velocity 
ballistic limit 
H. CONCLUSIONS OF BALLISTIC EFFECTS 

The multilevel ballistics analysis in this thesis shows that both design projectiles 
can successfully perforate the 100mm RHA target and the 200mm DLRC target. These 
results meet or exceed the ballistic design expectations of the railgun projectiles. The 
minimum expectation for the RHA target was penetration, and both projectiles clearly 
perforate it. An IFV projectile that is able to perforate 100mm of RHA at 2 km is at the 
forefront of other similar caliber armor piercing projectiles. Based upon the AUTODYN 
simulations for reinforced concrete impacts, the design projectiles will create perforation 
holes no smaller than their primary diameters. The primary diameters of the design 
projectiles are larger than the maximum primary diameters of American IFV armor 
piercing projectiles, which are no greater than 10 mm. 

I cannot conclude, based upon single shot perforation thicknesses, that the railgun 
projectile is more effective than an American IFV armor piercing projectile at creating 
the largest possible breach hole in concrete. The Bradley Fighting Vehicle has a 
sustained rate of fire of 100 rounds per minute; whereas, the railgun is only expected to 
achieve a rate of fire of 10 rounds per minute [11]. One would have to conduct a 
physical experiment between the two systems, allowing for sustained rates of fire, in 


order to decide which projectile had better overall perforation effect against concrete. In 
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conclusion, the proposed design projectiles demonstrate competitive ballistics effects 


against realistic Infantry battlefield targets. 
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Hi. AERODYNAMIC STABILITY AND DRAG 


A. AERODYNAMICS OVERVIEW 

During the development of any unguided kinetic energy projectile, researchers 
must thoroughly analyze in-flight stability and drag to ensure their design meets intended 
trajectories and impact velocities. In the development of the projectiles for this design, I 
performed an aerodynamic analysis as the second step in a recursive design loop that also 
considered ballistic effects and electromagnetic launch. My aerodynamic analysis is 
purely analytical, and relies principally on methods presented in Tactical Missile Design 
by Dr. Eugene Fleeman and one of his primary references, NACA Report 1307 [28], 
[29]. After ensuring that my design met the constraints of the aforementioned sources, I 
successfully calculated and proved both pitch and yaw neutral stability with a static 
margin of at least two percent. Continuing to use Dr. Fleeman’s analytical methods to 


calculate the Coefficient of Drag, C,, I concluded my analysis by developing an 


equation of velocity verses flight distance for each projectile. This final equation created 
a continuity link between the ballistic analysis and the electromagnetic launch by 
equating a particular launch velocity to an impact velocity. 

B. CONSTRAINTS AND ASSUMPTIONS 


The approach that aerodynamic experts use for the evaluation of subsonic, 
transonic, and supersonic flow is different for each flow regime. My design of the 
Infantry Railgun Projectile is constrained to supersonic flow and to trajectories that are 
strictly line-of-sight. The design’s lowest calculated velocity, at the limit of the furthest 
possible trajectory, is 828 m/s (See Table 9). Given a constant speed of sound at 343 m/s, 
this velocity equates to a Mach number of @=2.4. The highest possible velocity in the 
design is constrained to the launch velocity of 1100 m/s, resulting in a Mach number of 
M=3.2 (See Chapter HI, Section D). The geometric shape of the design projectile is 


approximately a right-circular cylinder with a right-circular nose cone. The length to 


diameter ratio for both projectiles is approximately 11.2 + 0.1 (See Appendix B2). 
Considering the geometric shape and the L/D ratio, the design projectile is a slender 


body, and is applicable to slender body theory [28]. My only aerodynamic wings in both 


af 


design projectiles are triangular tail fins. These fins are not banked, in a fixed position 
after launch, and in-line with the body axis [29]. One should consider these fins to be 


slender because their ratios of thickness, ¢, to mean aerodynamic chord, C are much 


“MAC > 


less than one, ¢/C\,,- U 1 (See Section C.3.) [28]. The aspect ratio, A, for these fins is 


approximately 3.85 + 0.02 (See Section C3). This would allows one to apply slender 
wing theory because the wing is thin and the aspect ratio is low or on the order of A=3. A 
problem with the application of slender wing theory is that it is only relevant for Mach 
numbers less than 1.53. The more appropriate theory is linear wing plus Newtonian 


impact theory. The constraint of this theory is satisfied by the following relationship: 


M > J1+(8/7A) ; where M is the Mach number and 4 is the aspect ratio [29]. 


Evaluating the right hand side of the argument, one should find a value of 1.2, which is 
lower than the minimum design Mach number of 2.4. From the discussion above, one 
can realistically constrain the aerodynamic analysis to slender body theory and linear 
wing plus Newtonian impact theory in the limit of the supersonic flow. 

Throughout the aerodynamic analysis, I had to make certain assumptions about 
the performance of the design projectile. For instance, while the projectile is in the 
electromagnetic launcher, its tail fins are folded inside the body. One must assume that 
these fins reliably deploy upon launch and lock into position. Although this design 
includes possible geometries and locations of fin-deployment pistons, I have intentionally 
emitted discussion of reliable deployment mechanisms. I also assume that the projectile 
does not roll upon launch or during flight. Although the current tail configuration is 
prone to an unstable rolling moment derivative, I would have to perform a wind-tunnel 
test or simulate the projectile in a finite element model to accurately analyze the rolling 
moment [28]. Taking a simplistic approach, I assume that the railgun does not impart 
spin to the projectile during launch, and that the projectile continues to travel to its target 
without any roll. One should consider these reasonable, but simplistic, assumptions 
given the vertical and horizontal plane symmetries of the railgun projectile. Since 
aerodynamic flow is symmetric in each of these planes, the appropriate stability fins can 
act in their respective planes to counter the nose-body normal force. The last assumption 


that one needs to make is that flow interference between the body cylinder and tail fins, 
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as well as between the yaw fins and the pitch fins, has no significant impact on stability 
calculations. In Tactical Missile Design, Dr. Fleeman states that for conceptual design 
analysis it is reasonable to isolate the normal force of distinct geometries on the projectile 
[28]. Although he agrees that interference calculations are needed for a more thorough 
analysis, he states that in conceptual design one can assume wing-body or wing-wing 
interferences are small compared to the effective normal force of the isolated geometries 
[28]. Overall, the assumptions that I discuss above are appropriate and necessary for the 


concept-level of analysis that I performed in this research. 


OF AERODYNAMIC STABILITY PARAMETERS 

In order to analyze pitch and yaw stability for the design projectile, I had to first 
calculate several aerodynamic parameters from the geometric specifications of my 
design. The theoretical development and specification of these parameters are found in 
Tactical Missile Design by Dr. Eugene Fleeman [28]. A secondary reference which Dr. 
Fleeman cites and that I also use is “NACA Report 1307” [29]. 

1. General Body Parameters 

a. Surface Reference Area 


The reference surface area, S,,-, of an aerodynamic body is the cross 
sectional area with respect to flow. For the case of a long rod projectile that has a 
cylindrical geometry with diameter D, the reference surface area is: Sp, =2D°/4. For 


cylindrical geometries, the reference surface area is also the reference body area 


(Sper =Spopy)- The design projectile falls into the category of a right circular cylinder, 


with the exception of its shaved sides. In Figure 19, one can find the geometric 
relationships needed to account for the exact cross section area of the design projectile. 
The areas denoted by Al and A2 are exactly equal and represent the missing cross- 
sectional area due to the projectile’s shaved sides. The parameter D1 is the primary body 
diameter of the projectile. Parameters W1 and W2 represent the body width and height 
with respect to the projectile’s shaved sides. I have specified the measurements of 
parameters D1, W1, and W2 for each projectile in the Appendix (See Appendix B2). 


Using the relationships for the reference area described under the Figure 19, I calculated 
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the following surface reference areas for the two projectiles: S',,,.(30mm) = 665.5mm* 


and S,.,(40mm) =1183.0mm?° . 


Design Projectile’s Body Cross-section 


ow 
(WW1) 42 


(W2) #2 


D1 = Projectile's Primary Diameter 

W1 = Height of Shaved Sides 

We = Thickness Between Shaved Sides 

A1=A2= Missing Area From a Circular Cylinder 

8 = Geometric Angle From Center to Shaved Sides 


Figure 19. Surface Reference Area 
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b. Equivalent Diameter 


One can calculate an equivalent diameter, ¢, for the design projectiles 
from the surface reference area. If one assumes the surface reference area represents a 
perfect circle, then the following relationship is true: S,,,=d°/4. Solving for the 
equivalent diameter, d, I calculated the following equivalent diameters for each 
projectile: d(/30mm) = 29.0 mm and d(40mm) = 38.8 mm. 

c. Reference Chord and Center of Gravity 

The reference chord length, c, is the mean geometric centerline that runs 
from the nose tip to the base of the projectile. The two projectiles have the following 


reference chord lengths: c(30mm) = 325mm and c(40mm) = 433mm (See Appendix B2). 
40 


The center of gravity for each projectile is X ce (30mm) =139.7mm and X cc(40mm) 
=186.2mm (See Appendix B7). 


pa Nose-Body Parameters 


The nose cone of a slender aerodynamic body is the primary geometry that 
separates the flow of air around the rest of the projectile. At supersonic speed, the shape 
of the nose can significantly impact the amount of energy lost to the creation of a primary 
Mach cone surrounding the body. The optimum aerodynamic nose shape is a perfect 


right circular cone that comes to a sharp point at the tip. The nose-fineness, /, , which is 
the ratio of the nose length to the base diameter, /,,/d , is a parameter that represents the 


quality of the nose. The general rating of common nose-fineness ratios are high, 
medium, and low. The value of each rating is 5, 2, and 0.5 respectively. For the Infantry 
Railgun Projectile, where I have scaled the 40mm projectile from the 30mm projectile, 
the nose-fineness is 3.44. This is between medium and high nose fineness and is a good 
value for limiting drag-loss to the primary Mach cone. Since the nose is the most 
significant geometry that separates flow, it is also where the body center of pressure is 
located in the limit of small angles of attack. This circumstance accounts for the position 
of the normal lifting force on the body. Therefore, the center of pressure for the nose is 
also the center of pressure for the body. I describe the semi-empirical relationships for 
calculating the nose-body center of pressure and normal force coefficient with respect to 
angle of attack in the following subparagraphs. 

a. Nose-Body Center of Pressure 

For slender bodies with right circular nose cones, the center of pressure for 


low angles of attack is at least one body-diameter from the nose tip [28]. NACA Report 


1307 defines a more exact relationship for the nose-body center of pressure, (X, led and 
is given by the following: (X,,), =Ly [1- Ve / (arly )| , where L,, is the nose length, 


ry 1s the radius of the nose base, and V,, is the volume of the nose [29]. For a perfect 


right circular nose cone, this relationship reduces to: 


(Xep)y = Ly 1-F (arity (arity) =0.67L,. 
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Using a nose length of 100.0 mm for the 30mm projectile and 133.3 mm for the 40mm 
projectile, I calculated the following nose-body center of pressure for each projectile: 
(Xp) y[30mm] = 67.0mm and (Xp) [40mm] = 89.3mm . 

db. Nose-Body Normal Force Coefficient 

The nose-body normal force coefficient that I use in stability analysis is 
actually called the coefficient derivative, (C,,,), , due to the angle of attack measured in 
radians. For slender bodies, this coefficient for the nose-body is approximately 2 per 
radian. As with the center of pressure calculation, NACA Report 1307 specifies a more 


exact relationship. The coefficient is related to the reference surface area and base 
diameter by the following: (Cy, )y = 2274 [Srp =2d°/2S pp, . Using this relationship, 


the equivalent diameter, and the surface reference area for each projectile, I calculated the 


following nose-body normal force coefficient derivatives: 
(Cy) [30mm] = 2.0[1/ radian] and (Cy,)y [40mm] = 2.0[1/ radian]. 
3. Tail Fin Parameters 


a. Tail Fin Planform Area 


The schematics in Figure 20 show relevant parameter specifications for 


the tail fins. The exact dimensions for each set of tail fins are shown in Appendix BS. 


Surface Planform ae St = St1+St2 St= Cr* bi2 


(Delta Confi —o | 
b/2 


|| b 
St1 = Fin 1 Surface Area St = Fin Planform 
Cr = Root Chord St2 = Fin 2 Surface Area Surface Area 
Ct = Tip Chord 
b = Wing Span = a 
Yep = Outboard Center of Pressure ‘ _ Tapel pave heal 
t = Wing Thickness A= Aspect Ratio = 2b/[(1+A)Cr] 
Cmac = Mean Aerodynamic Chord Yep = (b/2)/(3-A) 
XLE = Leading Edge of Cmac Cmac = (2/3)Cr(1+A+A*2)/(1+A) 


Figure 20. Tail Fin Planform Diagram 


42 


Two axial symmetrical tail fins represent one surface planform area. 
Both design projectiles have one yaw fin (YF) planform and one pitch fin (PF) planform. 
All fins on the projectile are delta configurations. This configuration is characterized by 
a right triangle, with no tip chord or taper ratio (A = Ct = 0). The following table 


summarizes the parameters identified in Figure 20. 


Variable i 30mm Projectile 40mm Projectile 
Yaw Fin Pitch Fin Yaw Fin Pitch Fin 
“ 26.4 22.4 35.3 30.0 
255 oh < 34.0 28.7 


51 68 57.4 
a 8.5 at lis 257 


eS SSeS 
CMAC 17.6 14.9 23.3 20.0 


eee 673.2 481.6 1200 861.0 
271.4 304.0 361.8 405.4 


Table 3. Tail Fin Planform Parameters 


b. Tail Fin Center of Pressure 


The center of pressure for any planform varies with the Mach number, 1/7. 
The Mach number is the ratio of the dynamic flow velocity to the speed of sound, 
M=Urc. The design projectile is intended to fly at low altitudes where the speed of sound 
is approximately 343 m/s. The dynamic flow velocity is no greater than the launch 
velocity of 1100 m/s. From design constraints in the electromagnetic launch, I set the 
launch velocity at 1100 m/s (See Chapter IV). Using a dynamic flow velocity of 1100 
m/s and speed of sound of 343 m/s, I calculated a Mach number for both projectiles of M 
= 312; 
Using linear wing theory plus Newtonian impact theory, Dr. Fleeman specifies an 
empirical formula to calculate a wing planform’s center of pressure with respect to its 


mean aerodynamic chord. One can examine this empirical formula in the following 


equation: (Xop)i. = Cyc (AB -0.67)/(248-1)+ X,,; where A is the aspect ratio, is 


M°*-1, Cyc is the fin’s mean aerodynamic chord, and X,,,is the leading edge of the 
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Ciuc With respect to the nose tip [28]. Table 4 shows my calculations for each planform 


on both projectiles. It is important for one to note that for all calculations, / is a constant. 
This follows from the definition of # and the fact that I have assumed the Mach number, 
M, is at a constant value of 3.2. This assumption is a conservative estimate because the 
aerodynamic efficiency of a wing or fin decreases with increasing Mach number. If one 
assumes the launch velocity Mach number for all calculations involving the fin’s center 


of pressure, then this will produce the lowest normal force coefficient for the fins. 


Variable | Units 30mm Projectile 40mm Projectile 
Yaw Fin Pitch Fin Yaw Fin Pitch Fin 
A ---- 3.86 3.84 3.85 3.83 
XLE mm 271.4 304.0 361.8 405.4 
(Xcp)TF mm 280.0 311.3 313.3 415.3 


Table 4. Tail Fin Center of Pressure 


C. Tail Fin Normal Force Coefficient 


The tail Fin Normal Force Coefficient, (Cy,,.)-, has the same quantitative 
definition as(C,,.),, Where it is a derivative due to angle of attack. For linear wing 
theory, Dr. Fleeman specifies the following relationship for calculating this coefficient: 


(C.J 3 4/ VM? -1 = 4/ per radian [28]. _ Using the highest possible Mach number 
as a constant, M@=3.2, I calculated a value of 1.33 per radian, which I use in all future 
stability calculations. As I stated earlier, when we use the highest Mach number then we 
are making a conservative estimate because a wing’s aerodynamic efficiency degrades 


for increasing dynamic flow velocities. 


D. PITCH STABILITY 


In the pitch stability analysis of the design projectile, | must sum two moments 
that act in the vertical plane and that rotate about the center of gravity. Using Figure 21, 
one can reference both moments from the center of gravity, which I have intersected with 


the pitch axis of rotation. 
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Figure 21. Pitch Stability Concept Sketch 


The nose-body moment is formed by the effective normal force, 


(Cya)v (Spopy/Seer), and the corresponding effective moment arm, 1,/d 
=[Xo¢ —(Xep)y]/d [28]. As one can understand from the sign convention in Figure 21, 
the nose-body moment acts in the positive direction, +(/,),,. I can further simplify the 
effective nose-body normal force by evaluating the ratio between the surface areas, 
Spopy / Sage. The body surface area is equal to the reference surface area. Therefore, the 
ratio is equal to one, and I can use (C\,,,), to represent the effective normal force. The 
pitch-fin moment is formed by the effective normal force, (Cy,,) pp (Spe /Spzp). and the 
corresponding moment arm, 7,,/d = [X¢g—(Xep) pr ]/d [28]. This moment acts in the 


negative direction to oppose the nose-body moment. I can apply the negative sign 
convention to the pitch-fin moment based upon its direction of rotation, or I can let the 
relationship between the center of gravity and center of pressure produce the negative 
sign analytically. One must note that I have referenced the center of gravity, Xcc, and the 
center of pressure, (Xcp)pr, from the tip of the nose where X=0. Therefore, the 
relationship of the center of gravity minus the center of pressure will yield a negative 
sign. One can further examine the summation relationship with the following equation, 
and the input parameters with calculated moments are quantitatively given in Figure 22 


and Table 5. 
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Figure 22. Pitch Stability Projectile Diagram 


Projectile | Projectile 
Derivative with respect to a 
Derivative with respect to a (W@=3.2 
Derivative with respect to a (/=2.4 


Pitch-Fin Planform Surface Area 481.6 861.0 
Cross-sectional Surface Reference 665.5 1183 
Area 


mm | Pitch-Fin Effective Moment Arm 


mm_| Body Effective Diameter 
I/rad | Sum of Pitching Moments due to 
le of Attack, a, at Mach 3.2 
yt M,), | \/tad | Sum of Pitching Moments due to -2.78 -2.83 
Angle of Attack, a, at Mach 2.4 
Table 5. Pitching Moment Input and Results 


The concept sketch in Figure 21 shows an exaggerated schematic of a stable 
projectile. The amount of rotation that a projectile undergoes in the vertical plane or 
about the pitch axis must be extremely small in order for it to maintain its intended 


trajectory. As the projectile travels against the direction of flow, the nose-body pitching 
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moment creates a torque about the center of gravity that is opposed by the pitch-fin 
moment. If the nose-body moment is greater than the pitch-fin moment, then the sum of 
the moments will be a positive value. In this situation, without active control, the angle 
of attack, a, continually increases and causes the projectile to spin off trajectory. If the 
sum of the two moments is zero then the projectile will have neutral static stability. 
Although neutral static stability is sufficient for the conceptual design of a rocket- 
powered missile with maneuverable wings and fins, it is inadequate for an unguided 
projectile. The pitch-fin moment must be larger than the nose-body moment in order for 
the projectile to reduce its angle of attack and turn back into the direction of flow. Ifthe 
pitch-fin moment is just slightly larger, the projectile will exhibit oscillatory motion 
about the direction of flow. This type of motion will not allow a projectile to maintain its 
intended trajectory, and precision will be poor. 

Aerodynamic experts rate a projectiles ability to maintain static stability by 
calculating its static margin (SM). The general definition of the pitch static margin, 


SM, is the difference between the pitch aerodynamic center, (X ,.),, and the center of 


gravity, X.,, divided by the projectile’s reference chord length, ¢ , which is given by 
SM, hes od aX, ay c [30]. Just as the center of gravity represents the point at 


which the gravitational force acts on an object, the pitch aerodynamic-center, (X ,.)p, 


represents the point at which the combined normal lift-force acts on an aircraft or 
projectile in the vertical plane. In order for a projectile to have static stability and a 
positive static margin, the aerodynamic center must be behind the center of gravity when 
one references these points with respect to the nose. For that case, the aerodynamic 
center will have a larger magnitude than the center of gravity. Before I can calculate the 
static margin, I must first calculate the aerodynamic center. 

I can calculate the pitch aerodynamic-center by refining the equation for the sum 


of the moments. By replacing the individual effective moment arms with the combined 


total moment arm, [x Be= ee Veale d,1can represent the sum of the moments with the 


aerodynamic center. When I include the aerodynamic center in the equation for the sum 


of the moments, I can solve for it explicitly. The pertinent equations are the following: 
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Given the final form of the equation for the pitch aerodynamic-center and the 
relationship for the static margin, I can conclude the analysis of pitch stability. The 
results of the static margin calculations will define each projectile’s pitch stability. For 
both projectiles, the static margin is approximately 2 percent, which is lower than the 
acceptable static margin of 10 percent for ballistic projectiles [31]. One must also 
understand that the static margins I calculated for the design projectiles are for the highest 
Mach number. This result provides two pieces of additional information. First, I cannot 
set the launch velocity any higher than 1100m/s. A higher launch velocity will decrease 
the effectiveness of the tail fins and result in a lower static margin. At a speed of 
approximately 1200m/s, both projectiles have a total pitching moment of roughly zero. 
The pitch stability analysis has set an implicit constraint in the design for the maximum 
launch velocity. Second, the static margin will improve throughout the projectile’s flight 
as drag naturally slows it to lower velocities. For instance, the 30mm projectile’s static 
margin will improve to 7.5 % when it is traveling at its impact velocity of 828 m/s. The 
input variables and the results of the aerodynamic center and static margin calculations 


are summarized in Table 6. 


Variable | Units Description 30mm 40mm 
Projectile | Projectile 
Center of Gravity 186.2 
Aerodynamic Center, M@/=3.2 195.6 
Aerodynamic Center, =2.4 21072 
Reference Chord 433 
2.2 
SM, % | Static Margin at Mach 2.4 Be 7.6 


Table 6. Pitch Static Margin Results 
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E. YAW STABILITY 

As one can observe from Figure 23, the concept sketch for yaw stability is almost 
an exact representation of the pitch stability sketch from the previous section. In this new 
orientation, I must analyze two moments in the horizontal plane or about the yaw axis of 
rotation. With the assumption that one can understand the previous discussion on pitch 


stability, I will reference that section for the yaw stability analysis. 


Pad 
Yaw Aoas 
of Rotation 


Pitch Axis 
of Rotation 


Figure 23. Yaw Stability Concept Sketch 


The effective nose-body normal force coefficient, (Cy,)y (Ssopy /Srer )» and the 
corresponding effective moment arm, ",/d =[X.4 —(X¢p)y]/d, are the same magnitude 


and dimensions as one finds in the pitch analysis [28]. Whether the projectile deviates 
from the direction of flow in the vertical or horizontal plane, it will experience the same 
nose-body normal force with respect to angle of attack. From the concept sketch in 


Figure 23, one can see that the nose-body moment is positive, according to the clockwise 
sign convention, +(M , ), . The moment that opposes the nose-body, in the 
counterclockwise direction, comes from the effective yaw-fin normal force coefficient, 
(Cyp)ye (Se /Spep ), and its corresponding moment arm, 1, /d =[X¢q—(X¢p)yp I/d [28]. 


The summation of the two moments and the development of critical-parameters follow 
the same method that I covered in the pitch stability section. The yaw moment 


summation relationship is given by following equation: 
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The input parameters and resulting moment calculations are given by Figure 24 and 


Table 7. 
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Figure 24. Yaw Stability Projectile Diagram 


Projectile | Projectile 
Deavative with respect to B 
Derivative with respect to B (M= a2 
Deniative with respect to B (2.4) 


| Pitch-Fin Planform Surface Area | Fin Planform Surface Area | 3673-2. | 2 1200 


ieee SS Surface Reference a 5 1183 
Rae 


| See | 
Nose-body Effective Moment Arm Hal 7 96.9 
ee Pitch-Fin Effective Moment Arm -140.3 -187.1 
a Body Effective Diameter 38.8 


- 1.51 


y 


- 3.90 


Table 7 Yaw Moment Inputs and Results 


During the pitch stability analysis, I developed and defined the relationship 


needed to calculate the pitch static-margin. I use the same relationship for evaluating the 
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yaw static margin, which is given by the following equation: SM, = te se pe |/e; 


ce | 


where SM, is the yaw static margin, (X,..), is the yaw aerodynamic center, X,,, is the 


projectile’s center of gravity, and Cis the projectile’s reference chord length. 
Accordingly, I need to calculate the yaw aerodynamic-center. Although an axial- 
symmetric missile or projectile will have the same aerodynamic center in different 
evaluation planes, this is not true for projectiles that have irregular geometries or non- 
symmetric planforms. The design projectile is not axial-symmetric and does not have 
symmetric planforms. In the horizontal evaluation plane, the total yaw moment is larger 
in magnitude than the total pitch moment in the vertical plane. The design projectile has 
a yaw static margin of 4 % at Mach 3.2. At Mach 2.4, the static margin increases to 9 %. 
As I stated in the pitch stability analysis, aerodynamic experts consider static margins 
above 10 percent as acceptable [31]. The relationship I used to calculate the yaw 
aerodynamic center is given by the following equation. The input parameters and 


resulting calculations for the yaw static margin are given in Table 8. 


dX(Mo), 


Ye 
“= keer + (Cys ee (S\;. Ses )| 


Projectile | Projectile 
X, mm _ | Center of Gravity [5907 186.2 
Aerodynamic Center, /=3.2 203:7 
Aerodynamic Center, /=2.4 225.5 
Reference Chord 433 
Static Margin at Mach 3.2 40% 
SM, % Static Margin at Mach 2.4 9.0 % 


Variable 


Description 


Table 8. Yaw Static Margin Results 


F, AERODYNAMIC DRAG 
The evaluation of the effects of drag on a projectile is necessary to determine 
realistic target impact velocities from an initial launch velocity. I used drag analysis 


techniques from Dr. Fleeman’s Tactical Missile Design to develop an equation for the 
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coefficient of drag with respect to velocity at a zero angle of attack [28]. Using this 
relationship, I then conducted a curve fitting analysis to represent the coefficient of drag 
as a simple function of a constant divided by the velocity. Finally, I used fundamental 
physics relationships to create an equation of each projectile’s velocity with respect to 
flight distance. These equations are my desired outcome for the aerodynamic drag 
analysis. They link the ballistics effects to the electromagnetic launch by producing an 
estimated impact velocity for targets at specific distances, given a fixed launch velocity. 
The zero angle-of-attack drag force, D,, on a slender body projectile is the 
product of the coefficient of drag, C,, , the dynamic pressure, q,,, and the cross-sectional 
surface reference area, S,,,. The expression for the drag force is given by 
Dy =Cp,9.Srer [28]. For the design projectiles, the total coefficient of drag is the 
summation of contributions from the nose-body, (Cp, )gopy, the yaw-fin planform, 
(Cp, yr» and the pitch-fin planform, (C;, )p,. The total coefficient of drag relationship is 
given by C,, =(Cp, sony + (Co, ye + (Co, pr [28]. I develop each of these contributing 


coefficients in following paragraphs. The dynamic pressure takes the same form that one 
can find in the Bernoulli Equation, where is it the product of one half the density of the 


fluid, p/2, and the square of the dynamic flow velocity, U7. The specific expression is 
given by 4, =(p/ 2: [32]. For all analytical calculations, I use the density of air at 
sea level or p=1.20 kg/m’. This assumption is reasonable because the intended 
trajectories of the design projectiles are for line-of-sight engagements at sea level. The 
dynamic flow velocity, U_,, is the relative velocity of the projectile, «, along the flight 


path. For example, the dynamic flow velocity is the muzzle velocity at launch, and it 
decreases during flight because of drag. From the discussion above, one can start to see 
that zero angle of attack drag force is proportional to the square of the projectile’s 


velocity, D,«<u°. By the end of the drag analysis, I will show an approximation 


expression of the drag force, where it is proportional to the projectile’s velocity, D, «<u. 
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The body drag coefficient, (Cp, ) sony » 18 the sum of three components that involve 
the primary Mach cone, (C), )sopy nave, the pressure drag from the base, (Cy, ) puse cousr > 
and the natural drag from skin friction, (C), ) sony reicnow- The form of the expression for 
the body drag coefficient is given by (Cy,)gony =(Cy,)eovrmae +(Cp, )awsecossr 
+(Cy )aopyrercrion [28]. The first coefficient, (Cp )sopy mare» Tepresents the drag loss 


from the projectile’s nose cone as it separates the dynamic flow around the rest of the 
body. For supersonic velocities or Mach numbers greater than one, a primary Mach cone 
forms and represents the flow discontinuity between the exterior dynamic flow and the 
interior flow around the projectile’s body. We can represent this coefficient for 


supersonic Mach numbers with the following relationship involving the nose fineness, 


Ly /d and the Mach number, M: (Cy )sopy ware = 3-6/[(Ly/d)(M -1)+3] [28]. The 


second coefficient, (Cp )pisecoasr > Tepresents the drag loss from base pressure for an 


unpowered projectile. The parallel body flow at the end of the projectile does not 
immediately close behind the flat base; rather, a low-pressure area forms causing 
rearward suction, which contributes to drag. The relationship for this coefficient at 
supersonic velocities is simply 0.25 divided by the Mach number, 
(Cp, \nase.coasr = 9-25/M [28]. The last coefficient, (Cp, ) sony meicron > characterizes the 
skin friction that we are most accustomed to thinking about when analyzing drag. Dr. 
Fleeman represents this coefficient with the following _ relationship: 
(Cp, sopy erection = 9-053(L/d)[M ((q,,L)}" ; where L is the total length of the projectile 
[28]. A problem will occur if someone tries to calculate values using mks units because 
Dr. Fleeman explicitly uses British units in his calculations. One needs to multiply this 
function by a factor of 1.7 in order to correctly calculate values in mks units [33]. The 


correct form of the function that I use in this analysis is the following: 
(Cp, sopy reicnoy = 9.0% L/d)[M (q,L))". The description of this last coefficient 


concludes the development of the body-drag coefficient. 


The drag analysis for the yaw-fin planform and the pitch-fin planform 1s the same. 


I will use the description for a tail fin to develop this analysis. One can equally apply this 


a, 


analysis to the yaw-fin and pitch-fin without changing any forms of the relationships. 
The coefficient of drag for a tail fin is the sum of contributions from wave drag, 


(Cp, ranwave> ANd skin friction, (Cp )pur reicriov- The complete form of the relationship 
is given by (Cp Jou =(Co, ran wave +(Co, ran ericrion [28]. The first coefficient, 
(Cp, runwave> Tepresent the supersonic drag loss from the tail fin’s creation of a 


subordinate shock wave inside the primary Mach cone. Dr. Fleeman states that this 
coefficient is based upon modified Newtonian theory in which it incorporates the 
calculation of the pressure across the fin’s normal shock as a function of Mach number 
[28]. The analytical relationship for this coefficient has a complicated form. One can 
find this relationship and physical parameter definitions in Figure 25. The last 


coefficient, (Cp, )rurrricrioy » tepresents the skin friction. It takes a similar form to what I 


introduced for the body friction coefficient, mentioned previously. Just as with the body 
friction coefficient, it is necessary to multiply by a factor of 1.7 in order to adjust 
calculations for mks units [33]. The relationship describing the tail friction coefficient is 


also shown in Figure 25. 
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Figure 25. Tail Fin Drag Coefficients 
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The only input parameters, which I have not discussed thus far, are two tail-fin 
angles, the number of wings, and the specific heat ratio specified in Figure 25. The 
leading edge angle, Atr, describes the amount of sweep for a planform. The angle is 
referenced from the body cylinder’s surface normal vector at the planform’s leading 
edge. All fins in this design are swept at approximately 45 degrees. The leading edge 
thickness angle, dLE, corrects the coefficient of drag functions for the sharpness of the 
fins. I have not specified this value anywhere else in the design because of the 
complications with drawing that type of geometry. I have assumed a leading edge 
thickness angle of 10 degrees for all fins. The number of wings, n,,, corrects the 
coefficient of drag functions for the number of planforms. In the case of the yaw-fins and 
pitch-fins there is only one planform for each set of fins. I have further specified the 


number of wings in Figure 25. The specific heat ratio, v, is the quotient of the specific 
heat capacity for air at constant pressure, c, verses constant volume, c, , and is given by 


vy =c,/c, . For the sea-level flow regime, this value is vy =1.4. 


Using a Maple computer worksheet, I was able to develop an equation for the 


coefficient of drag as a function of the projectile’s velocity, Cp (v7) =(C), )ropy 
+H(Cy yr +(Cy, Jor [34]. I created a continuous plot of this function for velocities 


between 343 m/s and 1100 m/s. Finally, I developed a best fit function in the form of a 
constant divided by velocity. Both projectiles had the same plot for their respective 
coefficient of drag functions. This is a reasonable expectation since both projectiles are 
equally scaled. One can find the Maple computer code for both projectiles in the 
Appendix (See Appendix E). The effectiveness of the curve fits for both projectiles are 


shown in Figure 26. 
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Figure 26. Curve Fit of the Coefficient of Drag for Both Design Projectiles 


From the results of the curve fit, I developed an equation for each projectile’s 
velocity as a function of flight distance. Using basic Newtonian physics, one can say that 


the sum of the forces in the x-direction, the direction of flight, is equal to the negative 
drag force, >) F, = M(du,/dt)=-D, =-Cy,¢,,Sper =—(/2)u,Cp, Spee. The result of 
the curve fit gives us a coefficient of drag that is approximately equal to: C, =570/u,. 


If I substitute this relationship into the force equation then I am left with the following 


differential equation: (du, /dt) =—(285pSp,;/M )u,. One can manipulate the left side of 


. . . . : 1 du 
the equation by separating the differential into the following form: = 


ia =U, aes . This last step should enable one to see that the velocity will cancel 
dx fit G8 


from both side of the differential equation. If I take the differential dy to the right side 
and integrate both sides of the equation over the limits of the flight path then I will have 
the final form of the equation for velocity as a function of distance. I already stated that 
this equation links ballistics to electromagnetic launch. The constrained launch velocity 


is 1100m/s, and from the ballistics analysis, I intend to impact targets at 2 km and 500 m 
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(See Chapter II). Given these input parameters and the velocity equation, I can solve 
explicitly for the impact velocity. One can examine the integration limits, the final form 
of the velocity equation, and a summary of results in the following equations and Table 9. 
In Table 9 I list two sets of impact velocities for each target. I previously calculated a 


coefficient of drag that had the form of C,, =525/u,. This value produced slightly 


higher impact velocities. For my ballistics analysis, I used these velocities as input 
parameters in my ballistic simulations. Considering the effectiveness of the Infantry 
Railgun Projectile against both targets, I consider the influence of this miscalculation to 


be negligible (See Chapter IT). 


‘ 2850S pep | 


du, =- eee fas > u(t) =1100m / 5 282 PSxee 


x [mks] 


Projectile | Projectile 
Surface Reference Area | 6.65 x107 | 1.18x10° 


|  M___|_kg_ | Total Mass of Projectile | 150 | 3.55 
[ui km) | vs | RHA at 200m _| RHA at 2000 m = a 390 
Corrected 
| u(500m) | m/s _|DLRCats50o0m_——|_—'1032,— | 1048 
(500m) | m/s | DLRC at500m 1024 1043 


(Corrected) 
Table 9. | Estimated Impact Velocities from Drag Analysis 


G. AERODYNAMIC CONCLUSIONS 

The concept-level analysis that I conducted for the design projectiles in the area 
of aerodynamic stability and drag satisfactorily achieved the overall design constraints 
(See Appendix A3). Using Dr. Fleeman’s analytical methods, I successfully proved 
neutral stability for both projectiles with static margins of at least 2 percent. Although 
the static margin results are not within an acceptable range for ballistic projectiles, I am 
confident that certain design modification will yield acceptable results. For instance, one 
could replace both the pitch and yaw stability fins with a cruciform lattice-fin design. 
This type of planform is appropriate for both the design constraint of being folded inside 
the railgun and the high supersonic Mach numbers that one can find in my design. I 


Sf, 


successfully achieved the second objective that I set out to complete. My secondary 
objective was to create an equation of velocity with respect to flight distance; my curve- 
fitting method results in reasonable impact velocities. In conclusion, the aerodynamic 


analysis of both design projectiles shows favorable characteristics for future research. 
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IV. ELECTROMAGNETIC LAUNCH (EML) 


A. EML OVERVIEW 

The results of the combined ballistic and aerodynamic analysis require a muzzle 
velocity of 1100 m/s in order to defeat the RHA target proposed in this thesis. I 
considered the required muzzle velocity for the development of both the 30mm and 
40mm railguns. From the design constraints in Appendix A, the lengths of both railgun 
barrels must be set at 3 meters. Barrel lengths shorter than 3 meters will require higher 
effective currents to attain the necessary launch velocity; higher current are undesirable. 
The rail separations in both railgun barrels are fixed based upon the thicknesses of their 
respective design projectiles. Given these specifications, there are a range of possible 
solutions for the heights and widths of the conducting rails. One can reference the 
dimensions of height, /, and width, w, in Figure 27. Matched with appropriate effective 


currents, each solution to the rail geometry can produce an 1100 m/s muzzle velocity. 
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Figure 27. Dimensions for Conducting Rails 


The design constraint of maximum temperature rise in the rails dominated this 
analysis and limited the number of solutions (See Appendix A9). The cross-sectional 
area of the conducting rails was the dominant factor in meeting the temperature 
requirement. Considering realistic guiding surfaces based upon each design projectile’s 
primary diameter, I selected the smallest possible values of / for both the 30mm and 
40mm rails. I then found the most favorable designs for both railguns by specifying w. 
Using fixed values for the rail geometry, I obtained appropriate effective currents in each 


railgun. These values were less than the design constraint of 2 MA (See Appendix A7). 


So 


Next, I made rough estimates on possible melting conditions for both projectiles. I used 
multiple finite element computer models to improve my rough estimates. The results of 
my analyses show that both design projectiles paired with their respective optimized 


railgun barrels will theoretically allow for a launch velocity of 1100 m/s without melting. 


B. THEORETICAL CONSIDERATIONS 

In ideal electromagnetic launch conditions, given by constant current and solid 
electrical contact, one can make reasonable estimates of the electromagnetic force 
induced by rectangular conducting rails. Over the past half-century, railgun researchers 
have developed theoretical relationships that allowed them to make reasonable estimates 
for the inductance gradient and the electromagnetic force [7]. These researchers have 
also developed thermodynamic relationships to estimate the electrical heating for all 
major components of a railgun [7], [8]. In realistic experimental conditions, there are 
several explicit and implicit factors that can significantly change the result of a launch 
with respect to the theoretical estimates. These factors include the exact pulse shape of 
the current, the loss of solid electrical contact between the projectile and the conducting 
rails, and unusual armature configurations. In order to address these issues adequately, 
one would either have to perform physical experimentation or computer modeling that 
included a transient current analysis with sliding contact. Realistic experiments and 
computer models with sliding contact are outside of the scope of this thesis research, and 
I have left them to future work on this design. Throughout my analysis of 
electromagnetic launch, I only consider theoretical relationships that are applicable to 


ideal conditions. 


The underlying theoretical principal that fully describes both the magnitude and 


direction of the electromagnetic force inside a railgun is the Lorentz force law. This 


fundamental physics principal allows one to calculate the magnetic force, F,, on a 


mag ? 


projectile by integrating the cross product of the volume current density, J, and the 


magnetic flux density, 6. It is represented in the following mathematical relationship: 


Ts, = fl (J x B)dt ; where dt is the volume integrand that depends on the coordinates 
V 


system. To evaluate this integral over a simple projectile geometry, one would either 
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have to use a finite element code or assume that the current and magnetic flux density are 
both uniform and unidirectional. A fact of electromagnetic launch is that the current does 
not distribute itself uniformly as it passes through a projectile. | Furthermore, the 
magnetic flux density through the projectile is not uniform. As the current approaches 
the conductive portion of the projectile, it bends along paths of lower electrical potential 
and resistance. Circular “Eddy Currents” form near the interface of the conducting rail 
and the first cross contact surface on the projectile [7]. Overall, the evaluation of the 
electromagnetic force on a projectile from the Lorentz force relationship is complicated. 
Railgun researchers have developed a single formula that is equivalent to the Lorentz 


force law, commonly known as the railgun force law [7]. 
fei (2 


/ is the effective current through the projectile. The inductance gradient, /’, has units of 
Henry per meter, or inductance per unit length. The value of the inductance gradient 
depends upon the geometry of the conducting rails, to include their separation distance. 
The Kerrisk method is a valid semi-empirical model for calculating Z’ in square or 
rectangular conducting rails [7], [35]. One can find the complete expression for the 


Kerrisk method in Figure 48 of Appendix C. 


Considering a constant or effective current, /, railgun researchers have used the 


railgun force law to develop a relationship for the conducting rail length, x, , in terms of 
an intended muzzle velocity, v, . The first step in obtaining this relationship is to express 
the force law in a form that describes the railgun projectile’s acceleration, a,. Since the 
mass of the projectile, M,, times its acceleration is equal to the railgun’s force, /, then 
a, =dv,/dt=L'I’/(2M,); where v, is the instantaneous velocity. If / is constant then 
the muzzle velocity, v,, is given by 


vy = fav, =< fla. 
iP: P 


0 
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Sincel=dQ/dt, then the charge that has passed through the projectile is 
Os Idt=2v,M,/(L'I). Solving for the projectile’s final position, x,, along the 
conducting rails at the muzzle velocity, v, , one finds 


L ' 
2M, 


iO 
4M, 


- a 
Nar JOM ge 


If one substitutes the expression for charge in terms of muzzle velocity, 


O=2v,M, [(L'l ), into the solution above, then the following final form for the rail 


length is given by 


This result contains parameters that can be used to solve for a required rail length. 
One can derive three subsequent relationships from the expression for required 


vi 


2 
| . If the rail length is a fixed parameter, then a simple 


rail length, x, - Hel 


manipulation of the relationship above will produce an expression for the final muzzle 


velocity, 


v, =LJ/L'x,/Mp . 


One can also rearrange these variables to produce an expression for the required effective 


current, 


T=v,JM,/(x,L'). 


Using the identity for current, / = dQ/dt , one can solve for the final launch time, /, . 


Q _ 2x, 


ee 
in| ta |aO= ; 


E 


All of these expressions, to include the Kerrisk Method, allow one roughly to match 


conducting rail geometry, projectiles, power supplies, and barrel lengths. 
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The average temperature rise in the conducting rails after an electromagnetic 
launch is a practical consideration when considering the ability of a railgun to sustain 
multiple shots. Railgun researchers have developed a heat relationship that allows for 
estimates on the average temperature rise, AZ’, in the conducting rails. This relationship 


is shown in the following expression: 
AT =| p./(2,C,)|G/A4’ ; 


where p, is the electrical resistivity, ¢ 


mM 


is the mass density, C,, is the heat capacity, 


G= } I’dt is the electrical action, and A is the cross-sectional area through which current 


flows [8]. The resistivity, density, and heat capacity depend upon the material properties 


of the conducting rails. The electrical action, G, is proportional to the electrical energy 


dissipated in a material in the form of heat, and has units of Ampere’ -seconds or 4’s | 


dv "7? 
From a BAe f , it follows that 
P dt 2M, 
v, =| dv, = is [Par= E C. 
J 2M, 2M, 
or 
gat 


7! 
This relationship for the electrical action, allows one to make reasonable estimates for the 
average temperature change in the conducting rails. 

Another practical consideration involving electrical heating is the ability of a 
projectile to withstand a maximum effective current over a particular launch time. There 
are limits on how much current can pass through conductive materials before they melt. 


One must consider two important parameters when evaluating potential melting situations 
in railgun projectiles: the electrical action, G = [7 “dt , and the cross-sectional area, A, of 
the solid conductive portion of a projectile that interfaces with the conducting rails. 


Using these parameters, railgun researches have developed a quantifiable material 


characteristic called the specific action to melt (SAM) [7]. One can calculate the specific 
action, g, of a conductive material by using the following relationship: g=G/A° [7]. 
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The SAM for aluminum, tungsten, and copper have respective values of 2.52x10*, 


2.43x10*, and 8.05x10* in units of | 4°s/mm* | [7]. Recent research on the analysis of 


electrical heating in conductive armatures has shown that aluminum and copper can be 
used right up to their SAM without diminishing the performance of their conduction 
ability or structural integrity [7]. The three materials mentioned above and their 
respective SAMs are important in this analysis because my design involves copper 
conducting rails and aluminum-tungsten projectiles. 

Each of the theoretical considerations mentioned above has specific relevance to 
the analysis that I conducted on the two design railguns. This analysis involves an ideal 


projectile launch with a constant effective current, /. 


OF BARREL OPTIMIZATION 

The temperature rise in the copper rails of my initial 30mm and 40mm railgun 
designs resulted in values near the melting point for copper. These initial designs 
involved optimizing the inductance gradient and the effective current. Successive 
attempts with this method consistently resulted in temperature changes above the design 
constraint of the 20 °C (See Appendix A9). These results indicated that the electrical 


heating in the copper rails was the dominate factor in this analysis. 


To eliminate any concerns about the temperature rise in the copper rails, I 


approached the design by setting the temperature change to 20 °C. The railgun 


relationship, AT = | Pp. if ( PC os G/ A’ , was used to give the temperature rise. One can 


eliminate the electrical action, G, in this relationship by substituting its ideal form, 


G=2M,yv, /L', to obtain 


Les eC, )| 2M,V, . 


AT = . 
aoe 


The ratio in the square brackets, consisting of the electrical resistivity, mass density and 
heat capacity, has a value of 0.005 [K—mm*/A’s] for copper. This value remains 


approximately constant for copper below its melting point [8]. The mass of the 
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projectile, MM, and the muzzle velocity, v,, are constant in both railgun designs. 


Applying the Kerrisk method to solve for the inductance gradient, L’, one will find that it 
is a function of the s, 4, and w (See Figure 27 and Appendix C). The values of s are 
constant in 30mm and 40mm railguns because they depend upon the respective thickness 
of their design projectiles. The last term in the thermodynamic relationship is the cross- 


sectional area, A = hw. 


During any optimization process, it is important to limit the number of variables. 
From the expression for A7 and the dependence of L’ on / and w in the Kerrisk method, I 
found that an increase in / past the minimum required contact surface for each projectile 
had a negative impact on the design solution. Thus w was selected to reduce AZ. In 
order to eliminate the rail height as a variable, I considered a practical launch situation 
involving the generic railgun barrel in this thesis (See Figure 27). Rails in the generic 
railgun are separated by ceramic insulators that also serve as a guiding surface for the 
projectile during launch. I made an assumption that a ceramic insulator extension of 2.5 
mm was sufficient to hold the conducting rails in position. This is a reasonable 
assumption, considering that the Lorentz force acts to push the conducting rails apart 
during a launch. I did not consider whether or not the 2.5mm ceramic insulator extension 
was sufficient to allow for tension on the rails before or after a launch. Given an overall 
extension of 5 mm, from above and below the rails, the minimum required / is the 
primary body diameter of the projectiles plus the extension. This results in 4 = 35 mm 
for the 30mm railgun and / = 45 mm for the 40mm railgun. With a fixed /, the 


optimization process was reduced to a single variable, w. 


I found the appropriate values of w for each railgun by adjusting this parameter in 
the thermodynamic relationship until the temperature change for the copper rails was 
approximately 20 °C. The resulting w for the 30mm railgun was 40 mm, and for the 
40mm railgun it was 50 mm. During the thermodynamic optimization, it was necessary 
to calculate Z’ as an intermediate step. Given the respective parameters in Figure 28, the 


Kerrisk method yields 2’ = 0.38 wH/m for the 30mm railgun and L’ = 0.39 wH/m for the 


40mm railgun. One can then calculate 7 and G from /=v,,/M,/(x,L') and 


G=2M,v,/L'. Having known values for L’ and /, one can estimate the resulting 
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Lorentz force on the projectiles using the railgun force law, F = L'I°/2. This concludes 
the optimization of the conducting rails. Figure 28 shows the configuration of both 


railgun barrels, and Table 10 summarizes the input values and the results of this analysis. 


30mm Railgun 


40mm Railgun 


Figure 28. Optimized Railgun Barrels 


Projectile | Projectile 
Launch time 
Him Inductance Gradient 


Variable 


Units 


Description 


7 


AT °C __| Temperature Change in Rails 
af MA | Effective Current 1.28 LOZ 
G A’s | Electrical Action 
F kN __| Force on Projectile 


Table 10. Inputs and Results for Railgun Optimization 


D. PROJECTILE HEATING 
With the temperature rise in the conducting rails constrained to a nominal level, I 
focused my analysis on the susceptibility of melting in the design projectiles. To 


estimate the effects of electrical heating, I compared my estimates for the specific action, 
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g=G/A’, of each material to its SAM. The railgun barrel optimization produced values 


for G, which are listed in Table 10. The cross-sectional areas, A, are difficult to define 
for these projectiles because of their complicated volumetric structures, as indicated in 
Figure 29. In order to apply the specific action relationship to this design, I had to make 
conservative estimates on cross-sectional areas. Furthermore, I had to analyze the 
conductive properties of the tungsten (W) and aluminum (4/) materials in order to 
appropriately adjust the electrical action to realistic values. My rough estimates indicate 
that g exceeds SAM for the aluminum armatures in both design projectiles. If the 
electrical energy dissipated as heat in the armature is averaged over the whole volume of 
the aluminum then g exceeds SAM for the 30mm projectile but not for the 40mm 
projectile. The difference between these two techniques indicates a need for further 


refinement. 


Side View 


[WWwOr] 
[Wwe] 


Tungsten Aluminum 
<$_——_-+ 


0 


[30mm] 59 60 immny 
Top View (40mm) 0 667 70 


Tungsten 


Figure 29. Conductive Portions of the 30mm and 40mm Projectiles 


A conservative estimate for the conductive cross-sectional area, A, in each design 
projectile consists of the W and A/ rail interfaces that have solid continuous volumes 
between the conducting rails. These areas are shown in Figure 29 as Al and A2, 
respectively. The vertical lengths of the two areas are equal. Therefore, one can use the 
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horizontal lengths to make a comparison. The length of the aluminum area is one-sixth 
or approximately 16 % of the total length along Al and A2. If W and A/ were to have 
similar conductive properties, then one might expect that approximately 16 % of the total 
current would flux through area A2, but in fact their resistivities differ. 

A comparison of the conductive properties between W and A/ is necessary to 
estimate the current distribution across areas Al and A2. The electrical resistivity, p,, of 
W is 5.7 pO-cm and for A/ it is 2.82 pO-cm [8]. At approximately 25 °C, 
p,(Al)<0.59,(W). Since the resistivity in a conductive material is a function of 
temperature, it is necessary to examine these values at the first melting point in the 
projectile. The melting point of A/ is 660 °C and the melting point of W is 3410 °C [7]. 
At the melting point of A/, the resistivity of W increases to 22.0 pO-cm and for A/ the 
resistivity increases to 11.2 wOQ-cm [8], [35]. Therefore, the A/ is still 50 percent less 
resistive or two times more conductive than W. 

Considering the area ratio and the difference in conductivity between areas Al 
and A2, one can make an assumption about the distribution of current. Since A/ is 
consistently two times more conductive than W up to its melting point, the distribution of 
current in A2 should be 32 % or 0.32 of the total current (i.e., /,/7 =0.32). Therefore, 
the remaining 68 % of the total current flows through the W in area A] (i.e., J, // = 0.68). 


The initial relationship that describes electrical action is the following: G= } I’dt. In 


this design analysis, the effective current, /, is a constant; therefore, one can assume that 
G«I°. This assumption allows for the development of effective electrical action values 
for areas Al and A2. These values are based upon the total electrical action in each 


railgun design and the square of the current distributions through each area. Therefore, 
the effective action in area Al is G,=G(J,/T) . Similarly, the effective action in area 42 
is-G, =G6(L/1 y. Using effective values for the electrical action and the specific areas 


of Al and A2, one can calculate the specific action for W and A/. A summary of the 


inputs and results for these calculations are shown in Table 11. 
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Variable Units Description 30mm 40mm 
Projectile | Projectile 


As Calculated Electrical Action 8.68x10° | 2.10x10"° 


Cross-sectional Area for W 1334 
Cross-sectional Area for A/ 266 
Effective Action for Al (W) 971x10° 


ts 5 Effective Action for A2 (A/) 8. ORIG! 2.15x10° 


Specific Action for Al (W) 11800 
sua Specific Action to Melt (W) 24270 
g2 2 + | Specific Action for A2 (A/) 29679 
SAM(AL) A’s/mm* Specific Action to Melt (4/) 25240 
Table 11. Inputs and Results for Projectile Heating 


The resulting specific actions for the A/ armatures, in both design projectiles, are 
higher than SAM(A/). The A/ armature in the 30mm projectile has a specific action that 
is 57 % higher than SAM(4/) , and for the 40mm projectile it is 18 % higher. This result 


indicates that both armatures will melt during launch. 


The conservative estimates that I made for the aluminum’s cross-sectional area, 
do not account for the entire volumes of the two armatures. Taking the area A2 for each 
armature and multiplying by their respective rail separations produces volumes of 3.9 
cm> for the 30mm projectile and 9.2 cm* for the 40mm projectile. In comparison with 


their total volumes in Appendix B3, these values account for 40 % of the armatures. 


Another approach to this analysis is to create effective areas for each material that 
are based upon their actual volumes. For the A/ armature, this method involves taking 
each armature’s total volume in Appendix B3 and dividing by their respective rail 
separation. This will an produce effective area of 381 mm for the 30mm armature and 
677 mm° for the 40mm armature. Using these values for cross-sectional areas in 
conjunction with their effective actions will produce specific actions of 27624 A’s/mm* 
and 21171 4°s/mm‘* for the respective 30mm and 40mm aluminum armatures. In this 
case, only the 30mm projectile shows indications of melting based upon a SAM(A/) of 
25240 A’s/mm*. I must assume that the experimental values for the armatures are 


somewhere between the original conservative estimates and this other approach. 
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One cannot theoretically determine the actual current distributions along the 
conductive portions of the design projectiles. The complex 3D geometries of the 
projectiles require finite element modeling in order to properly predict current 


distribution and electrical heating. 


E. EML FINITE ELEMENT MODELING 

I used Comsol Multiphysics, Version 3.2, to develop three finite element 
computer models for each design projectile [39]. The primary objective of this modeling 
was to simulate the current distribution and the time dependent electrical heating in the 
aluminum armatures. The projectiles are stationary in all models. A direct current (DC) 
model was used to establish the distribution of the volume current density. The solution 
to this model served as an input to a magnetostatics solver. I used the magnetostatics 
solutions to calculate the Lorentz force on the projectile. The output from the DC model 
was used again to obtain a final transient conductive heat transfer model and to estimate 
the incremental temperature change of each projectile during a 5 ms current flux. The 
results of the final heat analysis indicate that neither armature will melt under launch 


conditions. 


Computations of the electromagnetic launch included a portion of the copper 
conducting rails and an approximate version of the projectile for each design. The 
appropriate material properties of each component were used throughout the simulation 
process. I made the conducting rails long enough to contain the projectile. I placed the 
rails such that the aluminum armature was significantly more than four calibers past the 
current source, so that the computed Lorentz force could be compared with F = L'/?/2. 
Limitations with Comsol’s finite element mesh required a slight modification to the back 
end of the armature. The U-shaped curvature was approximated by four surfaces. A 
modeling environment surrounded the partial railgun structure in order to allow for 
magnetic flux density calculations. The geometric setup of each model is shown in 


Figure 30. 
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Modeling 
Environment. 


[m] 


Boundary Condition | ‘| 
(Electrical Ground) -- bipekhe 


“0.4 
Figure 30. 3D Finite Element Model 

The final six models, three for each design, were refinements to several 
previously investigated test models. The meshing of these models induced approximately 
7x10° degrees of freedom, and the computer that I used for simulations had only a 
maximum random access memory (RAM) of 4 Gigabytes. During the execution of the 
initial DC test models, I observed that current arced between the conducting rails. 
Several attempts at eliminating this problem involved adding sufficient and non- 
disruptive electrical insulation to prevent arcing, but this problem persisted into the final 
set of models. The current flux boundary condition that was supposed to induce an 
appropriate effective current in the test models did not create the same flux value at the 
projectile due to electrical arcing. However, I did observe a continuity of flux just after 
influx boundary and right before the electrical ground. I also observed a continuity of 
flux between both rails and the conductive contact surfaces of the projectiles. Therefore, 
I increased the value of the influx boundary to give the proper current through the 


projectiles. A summary of each model’s setup parameters is given in Appendix F. 
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i Direct Current Model 

The DC models for both designs show that current distributes itself throughout the 
tungsten and aluminum portions of the projectile. The modeling software allowed for 
post-processing of the final solutions. I was able to independently evaluate current flux 
through the contact surfaces of the aluminum armatures and the tungsten penetrators. 
Using Comsol’s volume integration postprocessor, I was also able to evaluate the 
resistive energy lost to the aluminum armatures. The volume integration evaluated the 
power dissipated to the armature in units of Watts. In order to find the exact energy in 


Joules, I multiplied the volume integration results by the time of launch, 7, = 5.45 ms, and 


used basic thermodynamic relationships to estimate the average temperature rise in the 
armatures. The results of my analysis indicate an average temperature in the armature 


that is significantly below the melting point of aluminum. 


Using Comsol’s boundary integration postprocessor, I evaluated the current flux 
on each contact surface between the conducting rails and two design projectiles. The 
current through the 30mm projectile was 1.27 MA. Approximately 445 kA (35 %) of 
current passed through the aluminum armature and 825 kA (65 %) passed through the 
tungsten penetrator. My rough estimates for the current distribution were 32% and 68%, 
respectively. The conductive portion of the 40mm projectile received a total current flux 
of 1.82 MA. The current distribution was 582 kA (32%) for the armature and 1.23 MA 
(68 %) for the penetrator. The 40mm projectile agreed exactly with the rough estimates 
for current distribution. Figure 31 shows the current distribution for the 30mm projectile, 
and Figure 32 shows a similar distribution for the 40mm projectile. The current 
distribution is shown with an arrow plot. Each arrow is a proportion of the magnitude of 
the current at its location. These figures also show the current densities on one x-y plane 


that cuts through the center of each projectile at z=0. 
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Figure 31. Current Distribution and Density in the DC Model (30mm Projectile) 
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Figure 32. Current Distribution and Density in the DC Model (40mm Projectile) 
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I evaluated volume integrals with Comsol’s postprocessor for the electrical power 
dissipated in the armatures. The dissipated power in the armatures was 668 kW and | 
MW for the 30mm and 40mm DC models, respectively. Multiplying these results by the 
theoretically estimated launch time of 5.45 ms, I was able to estimate the amount of 
energy, Q, lost to the armatures in the form of heat. This estimate shows that 3750 J of 
heat went into the 30mm armature, and that 5450 J went into the 40mm armature. The 
thermodynamic relationship for calculating the heat energy absorbed by solids or liquids 


is given by, Q=c,mAT ; where c, is the specific heat capacity, m is the mass, and AT is 


the average change in temperature. The specific heat capacity for the aluminum used in 
this design is 0.896 J/g-deg (See Appendix Bl). The mass of the 30mm armature is 
26.7 g, and for the 40mm armature it is 63.5 g (See Appendix B3). Using these values 
with the calculations for heat energy, I found that the A7 in the 30mm armature was 157 
°C and it was 96 °C in the 40mm armature. These values for AZ’ are well below the 
melting point of aluminum. 

2s Magnetostatics Model 

The purpose of the magnetostatics model was to test the validity of the 
distribution of current in the DC model. It used the output solution from the DC model to 
create a magnetic flux density, based upon the magnitude and direction of the volume 
current density within computational uncertainty. The resulting magnetic flux densities 
and field lines are in Figure 33 for the 30mm model and Figure 34 for the 40mm model. 
If the current distribution is correct, then the magnetostatics model should create a 
magnetic flux density that can induce a theoretical Lorentz force that is calculated from 


F=L'I’/2, with L’ from the Kerrisk method. The overall results for both railgun 
designs show agreement with L’ to within 6%. 
A quantifiable measure of validity is being able to produce the inductance 


gradient from the model’s calculation for the Lorentz force. Using the railgun force law, 


one can solve explicitly for the inductance gradient with the following expression: 
L'=2F/I’. Comsol’s postprocessor allows for a volume integral of JxB over any 
component of the railgun representation. Using the specific vector components of 
JandB that create a force in the launch direction, I evaluated the volume integration 
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function over the entire projectile. The results gave a force of 326 kN on the 30mm 
projectile and a force of 681 kN on the 40mm projectile. Again, the currents from the 
DC models were 1.27 MA and 1.82 MA, respectively. These remain unchanged and 
verifiable in the magnetostatics models. With the aforementioned values for force and 
current, I calculated L’ = 0.39 wH/m for the 30mm railgun and L’ = 0.41 wH/m for the 
40mm railgun. These values indicate differences between the theory and the model of 


2.6% and 5.1% within computational uncertainty for the respective railguns. 
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3. Transient Conductive Heat Transfer 

The rough temperature estimates from DC model were not specific enough to 
show hot spots or areas of concern for the two armatures. I used the output solution from 
the DC model and an initial temperature of 30 °C to setup a transient model in order to 
calculate temperature increases over 5 ms. I also used a geometry isolation function 
within the Comsol modeling software to prevent any heat exchange between the rails and 
the projectile. This modeling analysis shows that a large majority of both armatures 


remain at temperatures well below the melting point for aluminum. 


The required modification to the rear of the armature, in which I had to replace 
the smooth curvature with four surfaces, resulted in hot spots that exceeded the melting 
point of aluminum. It is reasonable to expect hot spots in situations involving current 
flow around edges. Overall, these hot spots only occur in roughly one percent of each 
armature’s total volume. Taking into consideration the number of hot spots, their 
maximum temperature, and their location with respect to the surface modifications, I do 
not believe that they properly indicate melting in the armatures. The location of these hot 
spots, near the rear of the armatures, does produce a concern about the possible melting 
of the PBI insulator, which has a maximum service temperature of 500 °C. A design 
modification that could mitigate this concern would be to place a thin heat shield between 
the armature and the insulator. One can gain a better perspective on the hot spots and the 
general heating of the projectile by examining Figure 35 for the 30mm projectile and the 


Figure 36 for the 40mm projectile. 
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Figure 36. 
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F. EML CONCLUSIONS 

The theoretical estimates and the temperature refinements from the computer 
modeling indicate that both design projectiles can conceptually achieve a launch velocity 
of 1100 m/s. This analysis shows that this velocity is attainable while satisfying the 
design constraints. There are specific details about the practicality of this design that 
need further consideration before full development and experimentation. This involves 
sliding contact and the achievability of the 40mm design’s effective current. The sliding 
contact between a tungsten penetrator and the copper rails may induce quick barrel wear 
in this application. The tungsten penetrator does not have a mechanism to maintain 
contact. It must rely on barrel tension and contact with the aluminum armature to sustain 
a current flux. The estimated effective current required for the 40mm projectile is close to 
the design constraint of 2 MA. One needs to consider the achievability of this effective 
current with respect to the Army’s version of the power supply for the FCS. From a 
general perspective, the design concept for the electromagnetic launch of the two 


projectiles demonstrates the potential for a future system application. 
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V. DESIGN CONCLUSIONS 


The research reported in this thesis shows the feasibility of designing a projectile 
suitable for use in urban warfare scenarios which does not produce potentially fratricidal 
components upon launch. The Infantry Railgun Projectile design concept creates a safer 
weapon for infantry forces. The most significant consequence of this self conducting 
projectile design is to create a deployable fin solution for aerodynamic stability. This 
thesis shows that the deployable fin design limits the maximum attainable velocity to 
1100 m/s. The ballistics and electromagnetic launch analysis indicates that this velocity 
still allows for both projectiles to meet their design requirements. Only through future 
work and experimentation will all of these issues find more validity. Overall, the design 
presented is a rough concept, but as this analysis shows there are aspects of the projectiles 


that should interest Army weapons developers. 


I was able to show the potential for design improvements throughout this analysis 
by maintaining a 30mm baseline projectile with a scaled 40mm version. One can identify 
how scaling affected each area in this design. In most weapons applications, the Army 
requires smaller and lighter equipment that performs just as well as larger and heavier 
versions. Between the two designs, the 30mm projectile is a more favorable choice 
because it met all the constraints of the design. However, there are areas where the 
40mm projectile performed better. With the same launch velocity as the 30mm 
projectile, the 40mm projectile had less erosion and a higher residual velocity after 
penetrating the RHA target. It also suffered less aerodynamic drag because of its mass. 
The 40mm projectile showed less electrical heating during electromagnetic launch, even 
though it required a higher effective current. The benefits of the 30mm projectile are of a 
practical nature. Its effective current is 660 kA lower than the 40mm projectile, which 
requires less pulsing by the Army’s CPA power supply. The 30mm projectile’s size and 
effective current will allow for a thinner and lighter railgun barrel because those 
specifications require less containment material. If the Army chooses to conduct future 
development work on this design, the comparison of the two projectiles presented in this 


thesis should facilitate informed design improvements. 
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APPENDIX 


A. DESIGN CONSTRAINTS 


1. The projectile penetrates medium armored vehicles 


The projectile penetrates 100mm of Rolled Homogeneous Armor (RHA) with a 
normal incidence obliquity angle at a distance of 2 km. Using my own investigation of 
contemporary medium caliber munitions, I assume that a successful penetration of this 
target by an armor penetrator demonstrates that it is effective against medium armored 


vehicles [17]. 
2. The projectile is multipurpose 


The projectile penetrates 200mm of Double Layered Reinforce Concrete (DLRC) 
with a normal incidence obliquity angle at a distance of 500 m, and creates a perforation 
hole diameter larger than 15mm. An ARL report states that a 200mm DLRC target 
represents a typical urban target [19]. In order for the design projectile to satisfy the 
perforation requirement, it needs to have a body diameter of at least 15mm. I derived the 
specification for the 15 mm perforation thickness from approximate body diameters of 
Bradley’s M793 TP-T training round, which is shown in Figure 37. Bradley IFV gunners 


believe this round is the most effective against concrete [11]. 


M793 TP-T 
BODY AND 
CASE, ROTATING 


PRIMER = CARTRIDGE tracer BAND ASSEMBLY 


ALUMINUM 
WINDSCREEN 


Figure 37. M793 (After Ref. [11]) 


3. The projectile is aerodynamically stable to provide precise trajectories 
The projectile must have at least neutral pitch and yaw stability. This minimum 
requirement will allow for aerodynamic modifications to improve the static margin to 


acceptable values. Static margin values = 10 % will create precise trajectories [31]. 
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4. The projectile minimizes SDZs for forward dismounted infantry 


No parasitic mass falls away from the projectile after launch. A projectile that 
does not eject any launch package components during flight will minimize Surface 
Danger Zones for forward positioned dismounted infantry. One can examine the impact 
of having fall-away launch components by examining the Army’s restrictions for the 


Bradley’s APFSD round in Figure 38. 


iy z i 3 
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Figure 38. SDZ Restrictions for a Bradley APFSD Projectile (From Ref. [12]) 


5. The projectile must have smooth and flat sides for the railgun 


The physical design of the projectile must incorporates geometries to maintain 
electrical contact with rectangular conducting rails. The circular cylindrical geometry of 
the projectile has shaved sides to fit the rectangular configuration of the barrel. The 
aluminum armature has an inverted “U” shape cut into its base to allow the 


electromagnetic force to push the armature against the rails. 
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6. The projectile must allow for feasible auto-loading 


Using an ARL concept for the FCS, shown in Figure 39, I have assumed that the 
maximum length of the ammunition storage compartment is 600mm. The ARL concept 
includes an autoloader. Therefore, the maximum length of the projectile can be no longer 


than 600 mm [2]. 


Ammunition 
Storage 
Compartment 


Figure 39. ARL Concept for FCS (After Ref. [2]) 


7. The railgun’s maximum average current is less than 2 Mega Amps 


The primary power supply candidate for the Army’s version of the FCS has one 
pair of Compensated Pulsed Alternators, shown in Figure 40. The individual peak 
current for each alternator is 1.08 MA. One could set the pulsing of the two alternators to 


get an effective current of approximately 2 MA. 


Pair of Compensated Pulsed Alternators (CPA) 


Peak Current per CPA = 1.1 MA 
Figure 40. ARL Concept for CPA (After Ref. [2]). 
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8. The conducting rails maximum length is no longer than 3 meters 


The barrel length restriction for the FCS is approximately 3 meters. This 
estimate is based upon a barrel that is roughly one-third the maximum vehicle length for 
the FCS, which is eight meters [5]. This FCS version comes from Armor Magazine and 


is shown is Figure 41. 


Future Combat System (FCS) 
External Dimensions and Features 
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Figure 41. Armor Concept for FCS (After Ref. [6]) 


9. The conducting rails can sustain 10 continuous shots without melting 


A proposed rate of fire for the Navy’s railgun design is 12 rounds per minute [9]. 
A naval railgun has the ability to eject heat into the ocean and can use large cooling 
systems. The Army’s proposed design for the FCS does not include a cooling system and 
has a limited ability to eject heat [2]. A 20 degrees temperature change per shot in the 
copper conducting rails will keep the average temperature of the rails below 500 °C after 
20 shots. The melting point of the aluminum armature is 652 °C, and the maximum 
service temperature for the PBI insulator is 500 °C. If the conducting rails reach 
temperatures above 500 °C, there is a possibility that the PBI insulator will melt while the 
projectile is being prepared for launch. These design characteristics limit the rate of fire 


to 10 rounds per minute for two minutes. 
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B. PROJECTILES 


The purpose of this appendix is to describe the 30mm and scaled 40mm design 
projectiles. The component overview identifies each major component with respect to 
design titles and material properties [15]. In following sections, I describe the geometric 
parameters and dimensional properties for each component. I design both projectiles 
with a Computer Aided Design software called Rhinoceros. This software allows for 
specific length and volume calculations [14]. This appendix ends with my calculations 
for the total mass and the center of mass of each design projectile. 


1. Component Overview 


Side View 


Titanium Stability Fins 


Tungsten Penetrator a nee Yaw & Pitch 


Tungsten, W2 a, 

Nensiy 193 glee. 

Tensile Strength: 980 MPa PBI Insulator ——_, 

Resistivity. .65C-6 chm-cm + gescit Celanese Celazole, PB 

Heat Cap: 0.134 Jig-deg_ Density: 1.3 g/cc 

Melting Point: 3370 ceg Tensile Strength: 160MPa 
Resistivity. 1E1/ ohrr cm 


. Max Service Air Temp: 500 degC 
Top View (Shorl Perivds) 


cence oo Aluminum, Al606°-T6 
er Density: 2.7 gicc 
2 on ae Tersile Strength: 310 MPa 
Z Resistivty: 3.99E-6 ohm-cm 
mo Heat Cep: 0.896 J/g-deg 
se Sa Melting >oint 652 degC 


Tilantiuint Alloy, Ti 

Density: 4.7 g/cc 

Tensile Streigth: 115GPa 
Resistivity. 1.7E4 ohm-cm 
Heat Cap: 0.525 Jig-deg 
Melting Point: 1600 deg C 


Figure 42. Design Projectile Component Overview 
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Ds Tungsten Penetrator 


SIDE VIEW 
Wi W3 L5 
. a id L3 (Red) 
L1 (nose) L2 (Body) Yaw Fin 
Armature aes Piston Base 
ae itch Fin 
TOP VIEW Tungsten Fiting: Slots Piston Base 
Center of Mass L6 


Omm 0 [mm] 133. 172.5 200 348 393 433 


Figure 43. Tungsten Penetrator 


Variable | Units Description 30 mm 40 mm 
Penetrator | Penetrator 
30 40 
10 13.3 
20 26.7 
100 133.3 
50 66.7 
L3 mm Rod Length 175 233.3 
L4 mm Yaw Fin Slot Length 28 37.3 
24 32 
2 27 
15 20 
Ww2 mm Body Width 26 34.7 
W3 mm Fin Slot Width 2 2.7 
V_W2 Total Volume 69.76 165.4 
M_W2 g Total Mass 1346 3192 
CM W2 mm Center of Mass 129.5 L725 


Table 12. Penetrator Dimensions 
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3. Aluminum Armature 


SIDE VIEW 


1 1 
1 1 1 


FRONT VIEW 


“Tie 
Jomm = 450-160 185. [mm] 
tl Distance with respect 
40mm 200 213 247 [mm] | fo Nose tip at O mm 
Figure 44. Aluminum Armature 


Variable 


Armature 
| D4 | =mm _[ConnectorDiameter_ | 3 


3 
Armature Length 35 


40 mm 
Armature 


Connector Length S 


Armature Thickness 


Shaved Side Height 15 


Body Width 
Interior Minor Width 


Interior Major Width 


Total Volume 


Total Mass 


Table 13. Armature Dimensions 
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4. Insulator 


Side View 


D1 


Insulator 
Center of Mass Li2 


40mm 200 213 325.5 348 393 433 [mm] 


Figure 45. Insulator Configuration 


Variable Units Description 30 mm 40 mm 
Insulator | Insulator 


D Primary Diameter 
D Rod Diameter 
D mm Base Plate Diameter 


1 

2 

3 

5 

L5 mm Pitch Fin Slot Length 

Base Plate Thickness 
Insulator Length 
Insulator Body Length 
Yaw Piston Length 
Pitch Piston Length 


mm Body Width 

Ww4 mm Interior Minor Width 
WS5 mm Interior Major Width 
M 
CM 


V_INS cm Total Volume 


mm Center of Mass 


Table 14. Insulator Dimensions 
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5. Stability Fins 
Side View Top View 


Pitch Fin 
Center of Mass 


30mm = 279.7 311.0 [mm] Distance with respect 
++ To Nose tip at 0 mm 
4 3729 414.7 [mm] 


Figure 46. Stability Fin Configuration 


Fins 
Diameter Yaw Piston 4 
D6 mm Diameter Pitch Piston 4 


Yaw Fin Length a1. 
Pitch Fin Length 32 
Yaw Piston Head : 11.9 
a 
36 
Yaw Fin Separation 36 
| Yaw Fin Thickness 1.3 
Pitch Fin Width 30.7 
Pitch Fin Separation 30.7 
1 Pitch Fin Thickness 1.3 
Total Volume, Yaw 0.898 
Total Mass, Yaw f 8.4 
Total Volume, Pitch 0.655 
Total Mass, Pitch E 6.2 
Center of Mass, Yaw 3129 


CM_PF mm ___| Center of Mass, Pitch 414.7 


Table 15. Stability Fin Dimensions 
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6. Total Mass 
M, = My, + M 4 Zi M ws + My of M pp 
30 mm Projectile: | Total Mass = 1497 grams or 1.50 kg 


40 mm Projectile: Total Mass = 3549 grams or 3.55 kg 


7. Center of Gravity 


s 1 
Xce = ae re oi M4CM + M wsCM ns a M,,CM,, 1 Mp-CM 5, ) 


T 
30 mm Projectile: Center of Gravity = 139.7 mm 


40 mm Projectile: © Center of Gravity = 186.2 mm 
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G RAILGUN BARRELS 

The purpose of this appendix is to describe the railgun barrels for this thesis 
research. The general concept of the barrel comes from an Institute for Advanced 
Technologies presentation by Dr. Mark Crawford [8]. I have slightly modified the IAT 
generic design to include a ceramic insulator. The configuration of the barrel and 


material properties are shown in Figure 47 [15]. 


Infantry Railgun Barrel 


Steel Confinement Material 
o= 700 MPa 
(Yield Strength) 


Confinement 
Thickness 


Ceramic xL = 3 m barrel 


Insulator 


Copper, Cu 
Density: 8.96 g/cc 


Resistivity: 1.7E-6 ohm-cm 
Heat Cap: 0.385 J/g-deg 
Melting Point: 1083 deg C 


Figure 47. Generic Rectangular Railgun Barrel 


The railgun force law, a function of inductance gradient and effective current, 


allows one to calculate the theoretical force on a projectile. This equation is shown in the 
following expression: F =L'/ of a where LZ’ is the inductance gradient and / is the 


effective current. Railgun designers commonly use the Kerrisk Method to calculate the 
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inductance gradient for a single pair of square or rectangular conducting rails [35]. The 


semi-empirical relationship for the Kerrisk Method is shown in Figure 48. 


Cross-section of Conducting Rails Kerrisk's Parameters 
A:= 044061 

B:= {-) 0.0771 
a1:= 3.397143 

az.= (-) 0.06603 
Ww b1:= 1.007719 
Vacuum Inductance Gradient Estimates , b2:= 2.743651 
Kerrisk's Method, LANL reportLA-9092-MS, b3:= 0.022093 
NOV 1981: b4:= 0.263739 


pe A+Bin{1+al™ +225) in| o1re2S+es% eae” | 
h hh hoo” COR 


Figure 48. Kerrisk Method for Calculating Inductance Gradients 
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D. BALLISTIC ANALYSIS 


I; AUTODYN 5.0 Modeling Set-up for MOD-2 Impacting RHA 

I used the following AUTODYN model setup for the 30mm and 40mm MOD-2 
projectiles impacting into the 100mm RHA target [25]. A general understanding of the 
AUTODYN graphical user interface is necessary to comprehend the setup summary. All 
input values for MOD-2 come from Figure 4 in Chapter II. All input values for RHA 


come from Figure 5 in Chapter IT. 


a. Modeling Domain 
-Select the two-dimension (2D) modeling perspective. 


{The x-axis is oriented from left to right, and the y-axis is oriented from bottom to top 
with respect to graphical images on a normal computer monitor} 


b. Load Materials 
[1] Tungsten 
-Select tungsten from the material library. 
-Ensure that the equation of state is “Linear-Von Mises.” 
-Change material density to 19.3 g/cc. 
-Change Bulk Modulus to 3.1E8 kPa. 
-Select Failure Model and input 0.005 for the Plastic Strain. 


{The Plastic Strain is approximately the Yield Stress, 1.51E6 kPa, divided by the 
Bulk Modulus, 3.1E8 kPa.} 


[2] RHA 
-Select RHA from the material library. 
-Ensure that the equation of state is “Linear-Von Mises.” 
-Change the Bulk Modulus to 1.7E8 kPa. 
-Select the Erosion Model and input 2 for the Erosion Strain. 


{A Century Dynamics AUTODYN technical advisor recommended the value of 2 
for the Erosion Strain [39].} 


-Select Failure Module and input 0.008 for the Plastic Strain. 


{The Plastic Strain is approximately the Yield Stress, 1.32E6 kPa, divided by the 
Bulk Modulus, 1.70E8 kPa.} 


ec. Create Parts 
[1] Block | and 2 (MOD-2 Rod Cylinder and Body Cylinder) 
-Select the Block geometry and fill with Tungsten. 
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{One must select individual Blocks for the Rod Cylinder and Body Cylinder. 
One can identify these geometries from the MOD-2 reference figure. The placement of 
block with respect to the x,y origin depends upon the preference of the user. The user 
must ensure that the Rod Cylinder Block is to the left of the Body Cylinder Block, and 
that the inner boundaries of the two blocks make contact.} 


-Input dimensions for length, x-direction, and width, y-direction. These values 
are based upon half-plane 2D axial symmetry along the x-axis. 


{All values for length, x, are the actual lengths that one can find from the MOD-2 
reference figure. All values for width, y, are the radial distances from the major-body 
axis or half diameters in the MOD-?2 reference figure. } 


-Input zoning that is two times the values entered for length, Ax=2x, and width, 
Ay=2y. 


-Input initial velocities which are the estimated impact velocities of MOD-2. 
{One can find the estimated impact velocities in Table 9} 
[2] Rhombus (MOD-2 Nose) 
-Select the Rhombus geometry and fill with Tungsten. 
-Select the Quarter Rhombus representation. 


-Input dimensions for the width, x-direction, and height, y-direction. These 
values are based upon dimensions for a full rhombus. 


{All values for width, x, are two times the actual lengths of the noses that one can 
find from the MOD-2 reference figure. All values for height, y-direction, are full 
diameters of the noses that one can find from the MOD-2 reference figure. The user must 
position the Rhombus such that the center body coordinates aligns with the bottom right 
edge of the Body Cylinder Block.} 


-Input zoning for “Across” that are equal to the diameters of the Body Cylinders 
from the MOD-?2 reference figure. } 


-Input initial velocities which are the estimated impact velocities of MOD-2. 
{One can find the estimated impact velocities in Table 9} 
[3] Join All 


-Select Join All for the Rod Cylinder Block, the Body Cylinder Block, and the 
Nose Quarter Rhombus. 


[4] Block 3 (RHA Target) 
-Select the Block geometry and fill with RHA. 


{The placement of this block must ensure that a 1mm gap is between the tip of the 
nose, which is to the left, and the bottom left corner of this geometry. } 


-Input dimensions for length, x-direction, and width, y-direction. These values 
are based upon half-plane 2D axial symmetry along the x-axis. 
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{All values for length, x, are the actual lengths that one can find from the RHA 
reference figure. All values for width, y, are the radial distances from the major-body 
axis or half diameters in the RHA reference figure. } 


-Input zoning that is two times the values entered for length, Ax=2x, and width, 
Ay=2y. 


-Input an initial velocity of zero for this geometry. 
d. Boundary Conditions 
-Create a Transmit boundary. 


-Select “Boundary Line” and place in the x-direction along the top of the RHA 
geometry. 


e. Rotate Model 
-Select “Rotate-270 degree” and “Axial-Symmetric.” 


{This action will take the current 2D half-plane axial symmetry and make a three- 
quarter 3D figure. All AUTODYN computations are evaluated in the 2D half plane and 
then mapped to the 3D axial symmetric figure} 


f Start Simulation 


2. AUTODYN 5.0 Modeling Set-up for MOD-2 Impacting DLRC 

I used the following AUTODYN model setup for the 30mm and 40mm MOD-2 
projectiles impacting into the 200mm DLRC target. Several input parameters in this 
model are exactly the same as the input parameters for the AUTODYN RHA model. I 
reference the RHA model in cases where there is exact replication by stating, “No 
Change from RHA Model.” A general understanding of the AUTODYN graphical user 
interface is necessary to comprehend the setup summary. All input values for MOD-2 
come from Figure 4 in Chapter II. All input values for DLRC target come from Figure 7 
in Chapter IT. 


Modeling Domain [No Change from RHA Model} 
Load Materials 
[1] Tungsten [No Change from RHA Model] 
[2] Concrete-L 

-Select Concrete-L from the material library. 

-Ensure that the equation of state is “Porous.” 
[3] Steel 4340 


-Select Steel 4340 from the material library. 
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-Ensure that the equation of state is “Linear-Johnson Cook.” 
Cc. Create Parts 
[1] Block 1 and 2 (MOD-2 Rod Cylinder and Body Cylinder) 
[No Change from RHA Model] 
[2] Rhombus (MOD-2 Nose) [No Change from RHA Model] 
[3] Join All [No Change from RHA Model] 
[4] Block 3 (Concrete-L Target) 
-Select the Block geometry and fill with Concrete-L. 


{The placement of this block must ensure that a 1mm gap is between the tip of the 
nose, which is to the left, and the bottom left corner of this geometry. } 


-Input dimensions for length, x-direction, and width, y-direction. These values 
are based upon half-plane 2D axial symmetry along the x-axis. 


{All values for length, x, are the actual lengths that one can find from the DLRC 
reference figure. All values for width, y, are the radial distances from the major-body 
axis or half diameters in the DLRC reference figure. } 


-Input zoning that is two times the values entered for length, Ax=2x, and width, 
Ay=2y. 


-Input an initial velocity of zero for this geometry. 

-Select the Fill function. 

{This additional Fill is to add one of the two Steel 4340 reinforcement bars. } 
-Select Ellipse and fill with Steel 4340. 


-Input an origin for the ellipse that is 26mm from the left edge and 145mm from 
the bottom of Block 3. 


-Input 6.5mm for the ellipse’s semi-major and minor axis. 

-Select the Fill function. 

{This additional Fill is to add the other Steel 4340 reinforcement bar.} 
-Select Ellipse and fill with Steel 4340. 


-Input an origin for the ellipse that is 177mm from the left edge and 145mm from 
the bottom of Block 3. 


-Input 6.5mm for the ellipse’s semi-major and minor axis. 
d. Boundary Conditions 
-Create a Transmit boundary. 


-Select “Boundary Line” and place in the x-direction along the top of the 
Concrete-L geometry. 
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e. Rotate Model 
-Select “Rotate-270 degrees” and “Axial Symmetric.” 


{This action will take the current 2D half-plane axial symmetry and make a three- 
quarter 3D figure. All AUTODYN computation are evaluated in the 2D half plane and 
then mapped to the 3D axial symmetric figure} 


x Start Simulation 


J; Critical Analysis Calculations 

The critical analysis theoretical method is a variant of the Tate Model for sub- 
hydrodynamic rod penetration. In this section of the appendix, I develop the primary 
relationships that one needs to make calculations with the critical analysis method. The 
primary reference for this analytical technique is the course notes from PH4857 [26]. 
The Department of Physics at the Naval Postgraduate School currently teaches PH4857 


as a primary course requirement for the curriculum of Applied Physics in Weapons. 
a. Cavity Propagation Speed and Critical Velocity 


[1] The first step in this analytical technique is to determine if 
the material strength of the projectile is larger than the material strength of the target, YP 
> RT. One can use the ultimate yield stress of a material to approximate the material 
strength. The yield stress of tungsten is 1.51 GPa, and the yield stress of RHA is 1.32 
GPa [27]. Using these values for material strength, we can conclude that the material 


strength of tungsten, YP, is in fact greater than the material strength of RHA, RT. 


[2] Given the case where our comparison of material strength 
satisfies the relationship of YP > RT, we can calculate the cavity propagation speed, VT. 
The relationship for the cavity propagation speed is based upon the pressure balance 


equation in Figure 18 and is approximated by the following equation: 


Ve= ———!——— ; where Vs is the impact velocity of the 


projectile, u is the square root of the ratio between the density of the target material and 


the projectile material (w= /p,/p,), and V- is the first critical velocity 
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(Vie? (Y,-R,)/p,). The input values for the cavity propagation relationship, and the 


results are shown in Table 16. 


Variable Description 40 mm 
Projectile | Projectile 
Vs Impact Velocity 890 m/s 
YP Material Strength of Projectile L531 


RT Material Strength of Target 1.32 

pp 3 | Density of Projectile (Tungsten) 19300 

pt > | Density of Target (RHA) 7850 
Target-Projectile Density Ratio 0.64 


m/s | First Critical Veloci 220 
m/s Cavity Propagation Speed 522.5 561.7 


Table 16. Cavity Propagation Speed Inputs and Results 


b. Penetration Depth at Critical Velocity 


[1] In order to create an equation for the penetration 
depth as a function of velocity, one must first consider the penetration depth as a function 
of time. We can assume that the penetration depth is limited by the cavity propagation 


speed. Therefore, the following integral can adequately approximate the penetration 


depth: P(t) = [V (eae We can manipulate this relationship to represent the 
0 


penetration depth as a function of velocity by making the following substitution: 

V(0)d0'=-2 1 WV, (w)dv"s where /(v’) is the un-eroded penetrator length as a 
P 

function of velocity and dv’ is the incremental change in velocity. The fundamental 

concept that allows one to make this substitution comes from the assumption that the 

elastic-plastic interface within the projectile is decelerating the remaining un-eroded 


penetrator. This assumption is summarized in the following force balance relationship: 
dv, ‘ ‘ , 
—-Y,A, = PplAp— where Ap is the cross-sectional area of the penetrator, / is the un- 


eroded penetrator length behind the elastic-plastic interface within the projectile, and 
dvp/dt is the deceleration. After making the appropriate substitutions, one will have the 


following relationship for penetration depth as a _ function of velocity: 
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Vo 
PY = = } l(v'\V,(v')dv'. One should note that this integral is specific to the limits 


PY; 


of the impact velocity, Vs, and the first critical velocity, Ve = V-. 


[2] One can derive the un-eroded penetrator length as a 
function of velocity from the force balance relationship inside the projectile at the elastic- 


plastic interface. The relationship for the un-eroded penetrator length is expressed in the 


lta -¥)} where / is the original 
P 


following relationship:  /(v')=Lexp 
201+yu 


length of the penetrator, and v’ is the instantaneous velolcity . 


[3] The form of the relationship for the cavity 
propagation speed does not change from what I have previously introduced. As a 


function of instantaneous velocity, the relationship for the cavity propagation speed is 


eS , 
Oe cad) 224 | cemaeeres Sear 2 
V(v')= is Ae) 


l+y by 


We can better 


expressed by the following: 


represent this function by considering a second order approximation for the instantaneous 
velocity. Given that the instantaneous velocity steadily decreases from the impact 
velocity, Vs, to the first critical velocity, Vc, it is approximately equal to half the impact 


velocity, v’ ~ (1/2)Vs. The final form of the instantaneous cavity propagation speed is 


2 


7 baste (1-0) . 


. vif 
l+u l-yw lt 


given by the following expression: V;(v') = 


[4] The final result of the integration for the penetration depth 
as a function of velocity, within the limits of the impact velocity to the critical velocity, is 
expressed in the following relationship. The input values and resulting calculations for 


the penetration depth at the critical velocity are shown in Table 17. 


PVG) ci: l—exp eee? 
Loseoug W1+pY,°* Fay 
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Variable | Units Description 30 mm 40 mm 
Projectile | Projectile 

Cavity Propagation Parameter 1.07 1.06 

:) 


kg/m 
GPa 


m/s_ | Impact Velocity 828 890 
m/s_| First Critical Velocit 220 220 
P(V,.)/L | ---- | Penetration Depth at Vc 


Table 17. Penetration Depth at Critical Velocity 


c Rigid Rod Penetration 


[1] After a penetrator reaches its critical velocity during 
impact, the resistance of the target’s material strength influences its deceleration. We can 
assume that the force balance equation at the elastic-plastic interface in the target is 


d 
approximately represented in the following expression: Pole ~ —R,A,; where /c is 


the un-eroded length of the penetration at the critical velocity. Starting with a function of 
rigid penetration with respect to time, we can manipulate this into a function with respect 


to instantaneous velocity. The integral form of this function and its final result are shown 


0 0 
in the following expression: P,,,, = fra —Pee vidv'= ee ys . From the 
Vo T 


T? Ve 
previous discussion of un-eroded penetrator length as a function of instantaneous 


velocity, we can assume that it has the following form at the critical velocity: 


1, = Lexp Ber Ss | —V2)|. The final form of the rigid penetration depth and the 
2(1+ L)Y> Hu 


resulting calculations are shown in the following expression and Table 18. 


vrei A exp] Sete) 


L 2R 2(1+ WY, : 


T 
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Variable | Units Description 30 mm 40 mm 


Projectile | Projectile 


kel m Dany of Projects (Tungsten) 19300 19300 


220 
Material Strength of Target ees 

Target-Projectile Density Ratio : 0.64 

Impact Velocity 890 

YP GPa_| Material Strength of Projectile 1.51 Rel 
Dipl ---- | Rigid Penetration Depth 0.071 0.055 


Table 18. Rigid Penetration Depth 


d. Total Penetration Depth 
[1] 30mm Projectile: a ae oy = —].24 
[2] 40mm Projectile: ws = a, 4 a8 —1.30 
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AERODYNAMIC ANALYSIS 


i Maple Computer Worksheet: 30mm Projectile 


Maple Worksheet: Coefficient of Drag Calculation for 30mm Infantry Railgun 
Projectile. * Note [ units are in meters, kilograms, seconds] 


> u:=u: 

c:=343: 
rho:=1.20: 
q:=1/2*rho*u%2: 
M:=u/c 

> L:=325E-3 

L n:=100E-3 
d:=29E-3: 

f n:=L_n/d 
gam:=1.4: 


GLE :=10*3.1416/180: 
TLE :=45*3.1416/180: 
MLE :=M*cos (TLE): 
Cmac_YF:=17.6E-3: 
Cmac_PF:=14.9E-3: 
tmac:=1E-3: 

b YF:=2*25.5E-3: 

b PF:=2*21.5E-3: 
Sref:=665E-6: 

S YF:=673E-6: 
S_PF:=482E-6: 
nw:=1: 

gam1:=gam/ (gam-1) : 
gam2 :=1/ (gam-1): 


> CDo_ body wave:=3.6/((L_n/d) * (M-1) +3) : 
CDo_base:=0.25/M: 

CDo_body friction:=0.053* (L/d) * (M/ (q*L) )*0.2* (1.71): 
CDo_body:=CDo_body wave+CDo base+CDo body friction: 


> CDo_YF_wave:=nw* (2/ (gam*MLE%2) ) * 

(( (gam+1) *MLE*2/2) “gam1* ( (gam+1) / (2*gam*MLE*2- (gam- 
1)))*gam2-1)* 

sin (dLE) *2*cos (TLE) *tmac*b YF/Sref: 

CDo_YF_ friction:=nw* (0.0133/(q*Cmac_YF)*0.2)*(2*S YF/S 
ref)*(1.71): 
CDo_YF:=evalf(CDo YF _wave+CDo YF friction): 
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> CDo_PF_wave:=nw* (2/ (gam*MLE%2) ) * 

(( (gam+1) *MLE*2/2) “gam1* ( (gam+1) / (2*gam*MLE*2- (gam- 
1))) *gam2-1)* 

sin (dLE) *2*cos (TLE) *tmac*b PF/Sref: 

CDo_PF friction:=nw* (0.0133/(q*Cmac_PF)*0.2)*(2*S PF/S 
ref)*(1.71): 
CDo_PF:=evalf(CDo PF _wave+CDo PF friction): 


> CDo:=CDo_body+CDo_YF+CDo PF: 


> £1t1:=570/u: 
plot ([CDo,fiti] ,u=343..1100,CD) ; 
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Ds Maple Computer Worksheet: 40mm Projectile 


Maple Worksheet: Coefficient of Drag Calculation for the 40mm Infantry Railgun 
Projectile. * Note [ units are in meters, kilograms, seconds] 


> u:=u: 

c:=343: 
rho:=1.20: 
q:=1/2*rho*u%2: 


_n:=133.3E-3 
d:=38.8E-3: 

f n:=L_n/d 
gam:=1.4: 


GLE :=10*3.1416/180: 
TLE :=45*3.1416/180: 
MLE :=M*cos (TLE) : 
Cmac_YF:=23.3E-3: 
Cmac_PF:=20.0E-3: 
tmac:=1.33E-3: 

b YF:=2*34.0E-3: 

b PF:=2*28.7E-3: 
Sref:=1183E-6: 
S_YF:=1200E-6: 
S_PF:=860E-6: 
nw:=1: 

gami:=gam/ (gam-1) : 
gam2 :=1/ (gam-1): 


> CDo_body wave:=3.6/((L_n/d) * (M-1)+3) : 
CDo_base:=0.25/M: 

CDo_body friction:=0.053* (L/d) * (M/ (q*L))*0.2* (1.71): 
CDo_body:=CDo_body wave+CDo_ base+CDo body friction: 


> CDo_YF_wave:=nw* (2/ (gam*MLE%2) ) * 

(((gam+1) *MLE*2/2) *gam1* ( (gam+1) / (2*gam*MLE*2- (gam- 
1)))*gam2-1)* 

sin (dLE) *2*cos (TLE) *tmac*b YF/Sref: 

CDo_YF_ friction:=nw* (0.0133/(q*Cmac_YF)*0.2)*(2*S YF/Sref) * 
(1.71): 

CDo_YF:=CDo YF _wave+CDo YF _ friction: 


> CDo_PF_wave:=nw* (2/ (gam*MLE%2) ) * 
(((gam+1) *MLE*2/2) *gam1* ( (gam+1) / (2*gam*MLE*2- (gam- 
1)))*gam2-1)* 
sin (dLE) *2*cos (TLE) *tmac*b PF/Sref: 
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CDo_PF_ friction:=nw* (0.0133/(q*Cmac_PF)*0.2)*(2*S PF/Sref) * 
(1.71): 
CDo_PF:=CDo PF _wave+CDo PF friction: 


> CDo:=CDo_body+CDo_YF+CDo PF: 


> £1t1:=570/u: 
plot([CDo,fiti] ,u=343..1100,CD) ; 
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F. ELECTROMAGNETIC ANALYSIS 


The purpose of this appendix is to describe the model setups for the 


electromagnetic launch simulations. The Comsol model setup summaries require a basic 


understanding of the software in order to properly replicate these simulations. All units 


in these models are standard metric (mks). 


1. 


Comsol Multiphysics 30mm Model Summaries 


a. Model Set-up 
[1] Application modes and modules used in this model: 


-Conductive Media DC (Electromagnetics Module) 
-Magnetostatics (Electromagnetics Module) 
-Heat Transfer by Conduction 


[2] Boundary mode: 


«— Modeling Environment ———____________, 


Boundaries(1-5,81] °  ¥ 


: oundary 241, Current Flux 


> Railgun ; 
“Boundaries [6-80)-" 
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[3] Subdomain mode 


[Mm] 
[4] Mesh Parameters 


Optimize quality 
Mesh geometry to level Subdomai 


Subdomain 


Element growth rate 


LZ 


[5] Mesh Statistics 


Number of degrees of freedom} 66424 


Number of edge elements 1683 
Number of boundary elements | 26146 


Number of elements 24393 
Minimum element quality 0.0784 


b. Application mode type: Conductive Media DC 


[1] Application Mode Properties 


Property 
Default element type] Lagrange - Quadratic 
Input property Fixed current density 


Weak constraints 


[2] Boundary Settings 


[Normal currentdensity(In) [OO 
[Normal current density (In) ]O fle [A/m2] 


[3] Subdomain Settings 


Electrical conductivity (sigma) 5.998e7 (Copper) 
Temperature coefficient (alpha) 0.0039 | 17e-6 (Copper) 
0. 


Se7 


0039 


Temperature coefficient (alpha) |0.0039 
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d. 


c. Application mode type: Magnetostatics (Electromagnetics Module) 


[1] Application Mode Properties 


Value 


agnetic 


Reference frame 
Weak constraints Off 


[2] Boundary Settings 


[3] Subdomain Settings 


E 
External current density (Je)| {'‘Jx_dcem';'Jy_dem’;'Jz_dcm'} 


Application mode type: Heat Transfer by Conduction 


— 
N 


[1] Application Mode Properties 


Property 
Default element type} Lagrange - Quadratic 
Analysis type 


Weak constraints 


[2] Boundary Settings 


1-5, 81 6- 
Thermal insulation 


ie.) 
fon) 


114 


[3] Subdomain Settings 


Thermal conductivity (k) 400 (Copper) 
8 


1300 
385 (Copper) 

[Heat source (Q) | 
Subdomain [89 MH2]10 
1300 
100 
Q_dem]Q dem [Q_dem 
Subdomain initial value 


Comsol Multiphysics 40mm Model Summaries 
a. Model Set-up 
[1] Application modes and modules used in this model: 
-Conductive Media DC (Electromagnetic Module) 
-Magnetostatics (Electromagnetics Module) 


-Heat Transfer by Conduction 
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[2] Boundary mode: 


“ ; ae en «= Modeling Environment 2 gf 
‘ Boundaries [1-5, 81] : y 


Boundary 21, Current Flux is x 


Boundary 16, Electrical Grouind 


Railgun ; 
Boundaries 6-80)" 


[m] 


[3] Subdomain mode: 


2 


M2 


b [Mm] 
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[4] Mesh Parameters 


[5] Mesh Statistics 


Application mode type: Conductive Media DC 


[1] Application Mode Properties 
Property 


Default element type 


Input property 
Weak constraints 
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[2] Boundary Settings 


[Normal current density (Im}O OT 
[Normal current density (ImJO | 9e8 [A/mm] 


9, 11-12 
1.77e7 
Temperature coefficient (alpha)| 0.0039|0.0039 |23.4e-6 


c. Application mode type: Magnetostatics (Electromagnetics Modules) 


[1] Application Mode Properties 


Static 


[2] Boundary Settings 


1-5, 81 
Electrostatic type 
Magnetic type 


[3] Subdomain Settings 


-12 


External current density (Je)} {‘Jx_dcem';'Jy_dem';'Jz_dcm’} 
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d. Application mode type: Heat Transfer by Conduction 


[1] Application Mode Properties 


Default element type] Lagrange - Quadratic 
Analysis type 


Weak constraints 


[2] Boundary Settings 


1-5, 81 6- 


Thermal insulation} Continuity 


oO 
jan) 


[3] Subdomain Settings 
Thermal conductivity (k) 400 (Copper) 

Density (rho) 1300  |8700 (Copper) 

Heat capacity (C) 385 (Copper) 


Heat source(Q) [0 

Subdomain [891210 
1300 
Q dem 
1-12 
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